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ABSTRACT 


Despite a huge technological advancement, air pollution remains one of the major 
problems in today’s world. Atmosphere is polluted by different gaseous species such 
as NOx (NO & NO 2 ), SO 2 , CO, and volatile organic compoimds (VOCs) etc. VOCs 
(e.g. benzene, toluene, and dichloromethane) have recently received greater attention 
in the field of environmental control due to its primary as well as secondary harmful 
impacts on ground air quality and human health. Due to its high BET surface area, 
activated carbon fiber (ACF) has been considered in recent times to be one of the 
promising adsorbents for controlling VOCs at low concentration level. 

In this work, an experimental set-up was designed and developed to carry out 
dynamic adsorption/desorption study on the removal of toluene at low concentration 
levels (< 1 % volume concentration) on the ACF wrapped over perforated teflon made 
reactor. The extent of removal of toluene was studied under various operating 
conditions such as, bed temperature (35-100 ®C),^as inlet concentration (2000-10000 
ppm), gas flow rate (0.2-1.0 slpm) and weight of the adsorbent (2-10 gm). A 
mathematical model was developed to predict the VOC breakthrough characteristics 
under these operating conditions. The model incorporates the effects of gas film mass 
transfer coefficient, pore diffusion and the adsorption/desorption rates within the pore. 

The experimental results showed that a bed temperature of 50 °C was most 
favorable for capturing toluene vapor from nitrogen stream. In a typical adsorption 
experiments carried out on 5 g of ACF sample the breakthrough and total adsorption 
times were found to be 7 and 50 min, respectively. The results also showed that the 
adsorption method was effective if the VOC emission levels were low i.e. parts per 
millions (ppm). The activated carbon fiber repeatedly showed a good regeneration 
capability on DC electrical regeneration. The experimental results and the 
corresponding simulated results from the model developed in this study were 
compared and foimd to be in good agreement. 
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CHAPTER-1 

INTRODUCTION 


Despite a huge technological advancement, air pollution remains one of the major 
problems in today’s world. Unhealthy air means unhealthy people. Air pollution has 
an enormous impact on everyone’s health, happiness and even the economy. Air is 
polluted by different gaseous species such as NOx (NO & NO2), SO2, CO, suspended 
particulate matters (SPM) and volatile organic compounds (VOCs) etc. VOCs (e.g. 
benzene, toluene, dichloromethane, and trichloroethylene) have recently received 
greater attention in the field of environmental control. 

Why? Emissions of VOCs have primary as well as secondary harmful impacts on 
ground air quality and human health. Eye, nose and throat irritation, headache, 
damage to liver and central nervous system, all may occur due to prolonged exposure 
to VOCs. VOCs may also have carcinogenic effects. Some of these organics can 
cause cancer in animals; some are suspected or known to cause cancer in humans. 
The role of VOCs as a precursor in the formation of photochemical smog and a 
number of toxic byproducts is well established. In atmospheric NOx-Os- 
photochemical chemistry, formation of hydroxyl (O//*) radicals is believed to be the 
main cause of chain reactions leading to the formation of O3. Atmospheric 
hydrocarbons and CO further aggravate the chain reactions leading to the formation of 
extra OH • radicals, which in turn produce extra O3 and many secondary byproducts 
such as HNO3 (source of acid rain), Peroxy acetyl nitrate (PAN) etc. The main steps 
involved in a typical NOX-O3-CO photochemical reaction are as follow: 

NOj +hu ^ NO + O 
O + O2+ M Oj+M 

O^ + NO^ NO2 + O 

O^+ho^Oi'D) 

0('D) + //j0^20//« 

C0 + 0H»-^C02+H02* 

HO2 • +NO -> N 02 isourceofexcess 02 ) + OH • 

OH • +NO2 -> HNOjiBy - products) 



What are the sources? In general, we can divide the sources into two categories: (a) 
stationary and (b) mobile. Organic chemicals such as solvents, thinners, degreasers, 
lubricants, and also liquid fuels are used and produced in many chemical industries. 
During handling, storage, and distribution of these chemicals, a significant amount of 
VOCs is emitted into atmosphere. VOCs can also be emitted from heating ventilation 
and air conditioning (HVAC) system (Wolkoff, 1995; Fanger et al, 1998) and 
bacterial and fungal activities (Rivers et al, 1992; Batterman et al, 1991). The 
extensive use of petrochemicals based organics and bio-organics in building products 
results in elevated indoor concentrations of VOCs. It has been seen that aliphatic and 
olefinic aldehydes and aliphatic acids are emitted from a large number of flooring 
materials, e.g. rubber, carpets, linoleum and vinyl (Wolkoff, 1995). Few polymeric 
materials like poly vinyl chloride (PVC) may also emit various types of VOCs. 

How to control? The VOCs emission control is presently one of the major activities 
in environmental program. Various control measures are taken by many industrialized 
countries to limit VOCs emissions into atmosphere within permitted levels. The 
commonly emitted VOCs in the atmosphere and the respective permitted emissions 
levels are listed in Table 1. There are various type of VOCs control methods presently 
available, including condensation, adsorption, catalytic oxidation and thermal 
oxidation, etc. ( Zesis, B., 1997). The adsorption of VOCs onto porous adsorbents, 
such as activated carbons, has been suggested as an innovative treatment process. 
Adsorption has been increasingly utilized in bulk separation or purification processes. 
The methods (especially condensation), other than adsorption are effective when 
VOCs concentrations are at higher levels (> 1%). However adsorption has been found 
to be suitable and economical if VOCs concentrations are at low levels i.e. parts per 
million (ppm) or sub ppm. Condensation followed by adsorption is found to be more 
effective if concentration of VOCs vary over a wide range (Verma & Gupta, 2001). 
Beside this, regeneration is a problem in other processes, however, in the case of 
adsorption; thermal regeneration is relatively easier to carry out. 

The keyword of an adsorption process is a porous solid medium. A high surface 
area or high micropore volume can be achieved due to the porous structure of solid 
and thus high adsorptive capacity. Steric, equilibrium and kinetic mechanism are three 
distinct mechanisms for the adsorption process. In steric mechanism the dimension of 
the pores is such that, all but the larger molecules can enter into the pores. In 
equilibrium mechanism, stronger adsorbing species is preferentially removed by the 



solid. The different rates of diffusion of different species into the pore are the main 
parameters in kinetic mechanism. Thus, faster diffusing species is preferentially 
removed by the solid. The success or failure of the process depends on how the solids 
perform in both equilibria and kinetics. A solid exhibiting good adsorptive capacity as 
well as faster kinetics is supposed to be a good adsorbent. Therefore in order to be a 
good adsorbent, a solid must have (a) reasonably high surface area or micropore 
volume and (b) relatively larger pore network for the transport of molecules to the 
interior. 

There are various adsorbents commonly used in industries such as alumina, 
silica gel, GAC (granular activated carbon) and zeolites, etc. These days, activated 
carbon fiber (ACF) has been considered to be one of the promising adsorbents for 
controlling VOCs at low concentration level. Due to its high BET surface area and 
pore size distribution (micropore volume), the adsorption capacity of ACF is higher 
than any other adsorbents. Pressure drop in the bed is an important factor. If the 
adsorbents are in powder (powder activated carbon, PAC) or granular (GAC) form, 
pressure drop in the bed is high, which is not desirable in industries. Regeneration is 
off course another important factor. In the case of ACF, regeneration can be done at a 
temperature of 1 50°C & at a lower flow rate, whereas a high temperature (350-400°C) 
and gas flow rate are required for the regeneration of GAC. ACF in bundled form 
becomes lighter and even stronger than steel. Carbon fiber is highly resistant to 
abrasion and chemical corrosion, and provides good lubricity. Not only as an 
adsorbent, ACF can also be used as catalyst supports, electrical and electronic 
materials for its unique pore structures and useful levels of electrical conductivity. 
Though it is a most complex solid, but due to its extremely high surface area and 
micropore volume it has been preferred by most of the industries. ACFs are 
synthesized from various homogenous polymeric base materials (precursors), such as 
poly acrylonitrile (PAN), phenolic resins, cellulose derivatives, and coal tar pitch 
[Mochida, 2000]. The essential microstructure or ACFs is developed during 
carbonization followed by activation of precursor materials by partial gasification in 
steam and/or CO 2 . The development of the microstructure is influenced by many 
factors, including the types of activating agent and operating conditions for 
carbonization and activation. Carbonization is essentially a de-volatilizing step to 
create pores within the materials. The subsequent step, activation causes the 
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development of a pore size distribution (PSD) and specific surface area responsible 
for adsorption. 

Following the introductory part and objectives, a review of the work done so far on 
the removal of VOC by adsorption on ACF by different investigators is discussed in 
the second chapter. A detailed description of the theoretical analysis and mathematical 
model, to understand the mechanism of the process is discussed in the third chapter. A 
description of the experimental set up is given in the fourth chapter. The fifth chapter 
describes the experimental results, along with the validation of the model developed. 
The conclusion drawn fi-om the present work and scope of the future work have been 
presented in the sixth chapter 


1.1 Objectives 

The objectives of this thesis work are as follows: 

1) Design and set-up of an experimental test bench to study adsoiption and 
desorptipon characteristics for the removal of VOCs on ACF. 

2) Determine breakthrough characteristics of the materials under v^ing operating 
conditions such as bed temperature, concentrations of VOCs, gas flow rate and 
weight of ACF, on the removal of VOCs. 

3) Study the effects of different operating conditions (bed temperature and 
concentrations of VOCs) on the removal of VOCs by two different samples of 
ACFs, one synthesized from viscous rayon & other from PAN (procured from 
abroad). 

4) Determine the optimal adsorption temperature and heat of adsorption for removal 
ofVOCsonACF. 

5) Carry out theoretical analysis and development of a mathematical model for the 
adsorption of VOCs on ACF, to understand the adsoiption mechanism and to 
predict the breakthrough characteristics. 

6) Model verification with the experimental data. 



CHAPTER 2 
LITERATURE SURVEY 

This chapter presents a critical review of various works done on the adsorption of 
VOCs by ACFs, together with different theoretical models developed to understand 
the mechanism of the removal process. An attempt has been made to identify the 
areas, which require further study. A review of various studies on the VOC adsorption 
by ACF reveals that most of the works may be categorized mainly in three groups: (1) 
equilibrium study to explain the isotherm data, (2) investigation of the effects of 
various operating conditions on the adsorption capacity of ACF for VOCs, such as 
temperature, concentration, etc. and (3) study to determine the breakthrough 
characteristics of VOCs on ACF under dynamic condition. 

Rood et al. (1994) have studied experimentally and theoretically, the 
equilibrium adsorption capacities of ACFs for acetone or benzene in air. The 
adsorption isotherms were measured and corresponding model equations were 
proposed for the ACF samples having three different BET surface areas. The 
concentration levels for the adsorbates chosen ranged between 10 and 1000 ppmv. 
The adsorption isotherm model proposed by Dubinin and Freundlich were used to fit 
the isotherm data. Dubinin-Redushkevich (DR) equation is one of the many equations 
developed by Dubinin to describe equilibrium adsorption of gaseous adsorbates. The 
equation is as follows; 

\nW = \nWf, -{Xq! j3cf ; where W and xo are the volumes adsorbed per unit mass 

of adsorbent [cm^/g] and slit-pore half-width [nm] respectively A good agreement 
between the model prediction and experimental data using the DR equation was 
achieved for both acetone and benzene. It was concluded for both the adsorbates that 
ACFs with the smaller slit pore half widths (related to lower BET surface areas) were 
more effective adsorbents than ACFs with larger mean pore sizes (higher BET surface 
area) at low adsorbate concentration levels (10-200 ppmv). 

Do et al. (1998) have studied a pure-component equilibrium adsorption isotherm 
of methane, ethane, and propane at different temperatures. The binary adsorption of 
methane in two heavier hydrocarbons (ethane and propane) by activated carbons 
having different pore structures were also studied. Six models were applied to the 
single and binary experimental data, namely: ideal adsorption solution theory (lAST), 
isoteric heat as a function of loading model (IHFL), micropore-size distribution 
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(MPSD), energy distribution (ED), and extended Langmuir and Sips models. lAST 
(Myers and Prausnitz, 1965) is based on a sound thermodynamic framework and 
widely used as a standard for multicomponent equilibria. In IHFL (Do and Do, 1997) 
model the variation of isoteric heat with loading is regarded as the measure of the 
heterogeneity of the adsorbent. This general model was tested with the single 
component adsorption isotherm data of many adsorbates on different samples of 
activated carbon and zeolites with good success. The MPSD model is based on the 
assumption that the surface heterogeneity of activated carbon is induced by the size 
distribution of the slit-shaped micropores. Sips (Sips, 1948) isotherm model is a 
combination of the Langmuir and Freundlich equations. This three parameter equation 
has been widely used for fitting the isotherm data of different hydrocarbons on 
activated carbon. It was found that all six models worked with a reasonably good 
accuracy to describe the experimental data of the hydrocarbon on commercial 
activated carbon. 

Verma and Dasgupta (2002) have studied the removal of VOCs from a gaseous 
mixture by cryogenic condensation, adsorption and a process combining both. 
Mathematical models were developed to predict the extent of removal of a binary 
mixture of VOCs in air by these two methods. The model developed to predict the 
breakthrough characteristics of VOCs consisted of an unsteady state species balance 
in the bulk gas phase, within the pores and on the surface of the adsorbing material. 
The kinetic rate expressions for the adsorption and desorption were obtained from the 
Sips equation. The Sips equation was obtained by equating the adsorption and the 
desorption rate. It was found that the Sips isotherm was suitable in explaining the 
isothermal adsorption/desorption behavior of a binary mixture of di-methyl chloride 
(DMC) and toluene. The model developed, was applied to simulate the experimental 
breakthrough curves reported elsewhere for the adsorption of various binary mixtures 
of organic compounds by activated carbon (Yun et al. 1999). 

In the following section a review on the effects of different operating conditions 
on the adsorption capacity of ACF for the removal of VOCs has been presented. 

Cal et al. (1995) have studied the effects of relative humidity (RH) on the 
adsorption of water soluble (acetone) and water insoluble (benzene) VOCs on 
activated carbon cloth (ACC). The RH values from 0-90% and organic concentrations 
from 350-1000 ppmv were examined. It was found that the presence of water vapor in 
the gas stream did not have appreciable effects on the adsorption of 500 ppmv 
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benzene until about 65% RH, when a rapid decrease resulted in benzene adsorption 
capacity with increasing RH. It was concluded that at the higher RH capillary 
condensation of water vapor occurred within the ACC pores, making them 
unavailable for benzene adsorption. It was also found that increase in benzene 
concentration might have a significant effect on the amount of water vapor adsorbed. 
At 86% RH and 500 ppmv, 284 mg water per g of ACC was adsorbed, while at 86% 
RH and 1000 ppmv, only 165 mg/g of water was adsorbed. It was also concluded that 
lower the benzene concentration the more profound the effect of water vapor on the 
adsorption capacity of ACC. The presence of water vapor in the gas stream with 
acetone had little effect on the adsorption capacity even at RHs of 90%. This is most 
likely because acetone is infinitely soluble in water, whereas benzene is insoluble and 
hydrophobic. 

Suyadal et al. (2000) have investigated the dynamic adsorption of 
trichloroethylene (TCE) in a laboratory scale packed-bed adsorber (PBA) on activated 
carbon at constant pressure (101.3 kPa). The experiments were carried out at different 
temperatures (25.6 °C < T < 35.8 °C) and with the TCE feedstock concentrations 
levels between 6300 and 8000 ppm. It was observed that TCE breakthrough curve 
shifted towards the left (short breakthrough tirne) with an increase in feedstock 
concentration. This behavior indicated a much quicker saturation of the adsorbent at 
higher concentration level. A shift in TCE breakthrough curve towards the right was 
also observed with decrease in temperature, which was attributed to an increase in the 
amount of adsorbed TCE. However no adsorption isotherm was proposed to fit the 
adsorption data 

Kang et al. (2001) have studied the adsorption behavior of polar methyl- 
ethyl-ketone (MEK) and non-polar benzene vapors on the viscous rayon-based 
activated carbon fibers in order to investigate the dependence of adsorption on 
adsorbate polarity. The influence of adsorbate polarity on the adsorption capacity and 
dynamic adsorption and desorption characteristics was studied. For polar MEK, it was 
found that a larger specific area resulted in a greater adsorption capacity if the 
concentration was greater than 100 ppm. However, reverse was the case when the 
concentration was less than 100 ppm. A larger specific surface area resulted in a 
higher adsorption capacity for non-polar benzene over the entire experimental range. 
The results also show that the adsorption capacity of MEK was greater than that of 
benzene for ACF with lower surface area. On the contrary, the adsorption capacity of 
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MEK was lower than that of benzene for ACF with larger surface area. It was 
concluded that the effect of adsorbate polarity on the adsorption capacity was more 
significant in the case of adsorbents with a smaller surface area. The time to reach 
equilibrium was much longer (2400 min) for non-polar benzene, whereas the time for 
MEK was only 1200 min. It was concluded that there was considerable oxygen 
complexes on the ACF surface, which had a greater affinity for the polar adsorbates 
and hence hastened the adsorption process. The desorption behavior of ACF with 
higher surface area was also examined by the adsorption of benzene and MEK at~ 
1000 ppm to reach equilibrium and subsequent heating to 393 K. The adsorbed 
benzene was totally removed in 0.5 h. However, it took 50 min. for MEK with ~ 2 
wt% still remaining on the surface. This suggested that a portion of the MEK is 
strongly bonded on the surface. 

Mangun et al. (1997) have studied the effects of pore size and pore volume 
on the equilibrium adsorption capacity and the adsorption kinetics for a series of 
normal alkanes in phenolic based ACFs. For the normal alkane’s series, the lower 
surface area (smaller pore size) materials had higher adsorption capacity for low 
boiling point gases/vapors and for adsorbing contaminants at low concentration 
levels. For higher boiling point gases/vapors and for adsorbing contaminants at higher 
concentration levels, the adsorption capacity was higher for the higher surface area 
materials (larger pore size) due to larger pore volume. 

Chiang et al. (2001) have studied the change of physical-chemical 
properties of activated carbons by ozonation and its effects on the adsorption of VOCs 
exemplified by MEK and benzene. The surface physical characteristics analyzed were 
BET specific surface area, micropore area, micropore volume and pore diameter. The 
BET surface area of oxidized activated carbon [AC (O 3 )] increased by 8.7% over that 
of the activated carbon (AC). The results indicated that oxidation increased mainly the 
micropore surface area, which contributed to 74% of the total increase in the specific 
surface area. Upon oxidation, the micropore volume also increased by 7.5%, along 
with an insignificant (0.5%) increase in pore diameter. Physical as well as chemical 
changes were also observed during ozonation. The results indicated that the 
concentration of oxygen radical within the functional groups, which remains on the 
surface of the activated carbon, slightly increased by 0.29%, whereas the same for the 
hydrogen decreased by 0.09% upon ozone treatment. The overall functional groups (- 
OH, -COOH, -C=0 etc.) were increased by 12%. It was concluded that the 



adsorption rate by AC( 03 ) is greater than that by AC, which could be attributed to the 
change of the physical characteristics of the activated carbon. The MEK adsorption 
capacity of AC (O 3 ) was greater than that of AC at various MEK concentrations. This 
was concluded as the surface of AC(0^) is more polar than that of AC and as MEK is 
a polar adsorbate, the MEK-specific adsorption sites will be more on AC( 03 ) than 
AC. 

Kang et al. (2002) have studied the breakthrough characteristics of low 
concentration MEK and benzene vapors in beds packed with rayon-based ACF having 
different surface areas. A smooth breakthrough curve (sigmoidal) profile was 
observed for ACFs having relatively larger surface area. However, at lower feedstock 
concentrations, the concentration of VOCs in the effluent increased gradually with 
time in the case of ACF having lower surface area, suggesting a large pore diffusion 
resistance during adsorption. This may be due to the narrow pore entrance in ACF 
leading to a low adsorption rate. Therefore, the ACF having lower surface area may 
not possess a good dynamic breakthrough adsorption characteristic, even though it has 
a higher equilibrium adsorption capacity at low concentrations. 

Larson et al. (1995) have developed a new ACC adsorption system that was 
integrated with electrothermal desorption and cryogenic vapor recovery. A laboratory 
scale fixed bed reactor was used to adsorb trace amount of VOCs on ACC. Once the 
adsorption was completed, the ACC was regenerated by electrothermal desorption 
under a small flow rate of N 2 gas. The resulting VOC gas stream was then treated 
(condensed) cryogenically with liquid nitrogen, making the pure toxic liquid available 
for reuse 

The breakthrough tests were performed on ACC and then with the Calgon BPL 
activated carbon in packed bed to compare the adsorption dynamics of the two 
adsorbents. The adsorption capacity results for Calgon BPL and ACC were 
determined to be 0.307 cm^/g and 0.553 cm^/g respectively. The breakthrough time 
for ACC and Calgon BPL adsorbents were 91 .7 min and 47.2 min respectively. It was 
concluded that ACC have higher adsorption capacity due to higher specific surface 
area for a higher specific pore volume. Electrothemical desorption was also carried 
out on ACC, which provided efficient regeneration and faster desorption rates. 
Stability for in situ electrothermal regeneration, control of pressure drop and easy 
handling were shown to be the other advantages of ACC in comparison to granular 
activated carbon (GAC). 



Verma et al (2002) have carried out both experimental and theoretical studies on 
the removal of VOCs by cryogenic condensation and adsorption. The adsorption 
experiments were mainly designed for the purpose of breakthrough analysis and 
carried out on both the ACF and activated charcoal under varying operating 
conditions to determine the relative performance of these materials with that of ACF. 

It was found that ACF exhibits larger adsorption for VOC than granular activated 
carbon under identical operating conditions. It was also observed that the 
breakthrough time during adsorption decreased with the increase in the VOC 
concentration, whereas condensation rate of VOC decreased with the decrease in the 
effluent concentration level. Hence, it was concluded, that adsorption was effective if 
VOCs concentration in the gas were in ppm and sub-ppm levels, whereas cryogenic 
condensation was found to be effective in relatively higher VOCs concentration levels 
(>1%). The regeneration of ACF was also carried out by electrical (~ DC 50 V and 1 5 
A) heating. A temperature range of 120-150 “^C and regeneration time 45-60 min were 
typically required for the complete regeneration of ACF 

Subrenat et al. (2001) have carried out the regeneration study on Joule 
effect in the case of ACC loaded with ethanol. The measurements of electrical 
resistance according to the intrinsic characteristics of the fabrics and temperature at 
the surface of the material for various activated carbon cloths have also been done in 
this study. The regeneration was carried out on ACC loaded with VOCs at a 
temperature of 100 ‘^C and 150 *^C and under nitrogen desorption flow velocity of 2.4 
m h’' . It was observed from the results that the increase in regeneration temperature to 
150 °C enhances the amount of desorbed compound leading to a desorption yield of 
100%. The desorption duration was also short (around 30 min) in comparison with 
that for lower temperature (100 ”C). The results showed that the electrical resistivity 
decreased linearly in the studied temperature range. It was also found that the 
electrical resistivity increased with the specific surface area and the activation with 
CO 2 results in a higher conductivity of the activatecT carbon fabrics, compared to 
steam activation. 

Hence, from the literature survey it is evident that quite a few works has been 
carried out on the equilibrium study to explain the isotherm data, whereas only a 
limited number of studies have been carried out to determine the breakthrough 
characteristics of VOCs on ACF under dynamic condition. The focus of the present 
work is on the application of regenerative type of sorbent, ACF in control of VOCs, 


10 



which is also a potential adsorbent for several other gaseous pollutants, such as SO2, 
NOx etc. In this thesis, both experimental and theoretical studies have been carried out 
to determine the breakthrough characteristics of the ACF during adsorption of VOCs 
under varying operating conditions. In addition, an electrothermal regeneration has 
also been studied for the desorption of VOC and making the ACF for reuse. A 
detailed description of the theoretical analysis and mathematical model, to understand 
the mechanism of the process is discussed in the next chapter. 



CHAPTERS 

THEORETICAL ANALYSIS 

In the present chapter, a mathematical model is presented to predict the time- 
dependent (unsteady state) concentration profiles of the adsorbing species (VOCs) on 
a solid adsorbent (ACF) under isothermal conditions. One important point that may be 
noted is that, in reality the steady state conditions never exists in the bed of adsorbing 
materials during adsorption/desorption. Hence, a finite adsorption/desorption rate 
always exists in the bed. The steady state is achieved only when the bed reaches 
saturation levels. The theoretical analysis in this work takes into account the effects of 
pore diffusion, in addition to the film mass transfer resistance around the fibers on the 
performance of adsorbents. The principal aspect of the present model is that it is 
based on non-equilibrium approach and thus incorporates the individual kinetic rate 
expressions for adsorption and desorption rather than equilibrium or pseudo- 
stationary assumptions made in most of the models (Yang, 1997) developed for 
adsorption in a fixed bed. By equating the adsorption and desorption rates an 
appropriate equilibrium isotherm such as Sips isotherm is obtained. 

The various assumptions made in the theoretical analysis in developing mathematical 
model for adsorption of VOCs on ACF may be summarized as follows; 

(1) Considering the fact that the concentration of VOCs are usually in ppm levels and 
an unsteady condition exists throughout a typical adsorption/desorption cycle, the 
effects of exothermic heat on the rate of adsorption were assumed to be negligible 
(Verma and Gupta, 2002). 

(2) Negligible pressure drop exists in the bed (Rood et al. 1995). 

(3) There is a constant fluid velocity throughout the bed. 

The present model is based on three governing equations: (a) species balance of the 
adsorbing component in the bed, (b) species balance of the component inside the 
pores of the fiber, and (c) surface adsorption within the pores. 

Species Balance of the component in the perforated bed: 

The experiments were carried out on ACF in cloth form wrapped over a vertical 
reactor made of Teflon. The outer surface of the reactor was perforated with one end 
closed. The reactor was encapsulated in a SS tube with provisions for gas inlet and 
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outlet. The VOC laden gas fed to the reactor came out radially into the SS tube. 
Consider a gas flow in a perforated bed of ACF under isothermal condition, in which 
pressure drop is negligible and there is no variation in fluid velocity in the radial 
direction the bed. By making a species balance in the gaseous phase in the ACF bed 
of porosity ,8, the equation is obtained as follows: 


dt ^ m R dR 




(3.1) 


The terms on the left hand side of Eq. (3.1) are transient, convection (radial direction), 
dispersion and mass transfer rate from the bulk gas to the fiber surface, respectively 
Here, 


C(} = Interfiber concentration of the species, at any arbitrary location r in the bed. 
Vr= Gas velocity in the radial direction 
Dr = Radial dispersion coefficient in the perforated bed 
a = external surface area per unit volume of the fiber 
= 4.0/ Di- 


Dr = Diameter of the fiber 


km = average mass film transfer coefficient around the fiber. 
C®F= Copcentration at the outer surface of the fiber 


Mass Balance of the component inside the pores of the fiber: 

Assuming the shape of the pores to be spherical, the concentration of the species, Cp 
in the pore is obtained from the following equation: 


dt dr ^ ’ 


In the above equation the terms in the left hand side are transient, diffusion and 
adsorption terms, respectively. 

Here, 

Cp = Molar concentration of species at a location r inside the pores 
a = Intrafiber void fraction 


a® = Total adsorption surface area per unit volume of the fiber 
Nr = Molar diffusion flux of the component in the radial direction 
Cs = Molar surface concentration of adsorbed species inside the pores 



Adsorption/descirptioii within the poresi 

The rate of change of surface concentration of the adsorbed species 
detennined from the individual surface rates of adsoiption and desoiption 


dt 


IS 


c 

dt “ " 


c 

l-(^) 


c 


A" max J 


- k A~)" 


c 


(3.3) 


A max 


Here, ka and kd are adsorption and desorption rate coefficients. 

At equilibrium, as the change in the surface concentration is zero, i.e. — ^ = 0, the 

dt 

corresponding isotherm for Cs can be obtained from equation (3.3) by equating the 
adsorption and desorption rates as follows. 


C, = C,. 




(3.4) 


Here K= adsorption equilibrium constant = 


k 

K 


The above equation is known as the Sips isotherm {Do et al, 1998) for single 
component system. As shown in equation (3.4), the Sips isotherm has three 
parameters; Csmax, K, and N, which are usually obtained from the experimental 
adsorption data under equilibrium condition. It may be pointed out that alternate 
kinetic rate expressions obtained from other isotherms including the extended 
Langmuir and Temkin can also be used to determine the rate of change of surface 
concentrations. In the present study, the Sips isotherm was selected due to the fact that 
the Sips isotherm has been found to explain the equilibrium isotherm data for a 
number of VOCs e.g. benzene, toluene and xylene etc (Do et al, 1998). 

The governing equations (3. 1-3. 3) are the basis of the model developed in this 
study for predicting the adsorbents performance in the bed. These equations are 
coupled^artial differential equations, with time (t) and radial directions (both R and r) 
as independent variables. Hence, in the present work an approach is adopted to 
simplify the numerical computations and reduce the CPU time. In this approach radial 
concentration profiles within the pores are averaged and the average surface and gas 
phase concentrations within the pores of the fiber are determined (Yang, 1997; Verma, 
2002). The detailed computational steps are shown below: 



First, the particle volume average quantities (Yang, 1997) are calculated as follows: 

(3.5) 

0 

— 1 

Q = \c./^dr (3.6) 

0 

Where C ^ and are the volume average surface concentrations in the pores and on 
the surface, respectively and Rf is the radius of the fiber. 

Second, by integrating equation (3.2) with respect to r and using volume averaged 
quantities and noting that = 0 at r=0 (due to symmetricity), the following equation 
is obtained; 


dC, 


■ P ^ \r 

a - + — N. 


dt R 


f 


dt 


(3.7) 


Where, N/ is the diffusional flux at the fiber surface. 

At the fiber surface, the diffusing flux into or out of the pores, N® is given by 

v;=-4.(c-c;) (3.8) 

Finally, to further simplify the model, a parabolic concentration profile within the 
pores of the fiber is assumed, which is a good assumption when compared with the 
profiles obtained from rigorous approach (Yang, 1997). Thus we have the following 
assumed expression for Cp. 

C^ = E + Fr^ (3.9) 

Where the constant E and F can be evaluated by using surface condition: 

=C/, at r = R^ 

E = C; - FR/ 

By integrating equation (3.9) and again using the volume averaged quantities, F is 
computed as follows, 

(3.10) 

2Rf 

C/ is determined by equating mass flux across the external gas film to that at the 
mouth of the pores within the fibers: 





fdC„) 


2D,FR^ 


(3.11) 


From the above expression Cj, is obtained in terms of independent variables Cp and 
C(i is as follows, 

C ' = KC\.+mDJD,)C\ (3 12 ) 

(k^+lODJD^) 

Here, De is effective diffusivity inside the pores 

Hence, by the approximation mentioned above equation (3.7) can be rewritten as: 


— — (C„ - C/) + — — = 0 
dt Rn ' a dt 


(3.13) 


In the present model, the simplified governing equations (3.1), (3.3) and (3.13) 
completely describe adsoiption/desorption behavior of a single component VOC 
(toluene) in a perforated bed of ACF in a gas flow under isothermal condition. 


Initial and boundary conditions for the above equations are as follows: 


At t=0, Cg= Cgi C» = C, 




(3.14a-c) 


where, Cg; is the concentration of trace ainount of VOC which may be present in the 


AtR = R,, C<, = Cc,f. C, =C.^ , „„ , C, =C<y 

where, Cor is the concentration of VOC in the feed stream 


(3.15a-c) 


At R — R'5 


^ = 0,^ = 0, 


(3.16a-K;) 


Non-dimensionalisation of governing equations and initial and boundary 
conditions: 

The following dimensional variables are introduced to render the governing equation 
dimensionless: 

Bulk Gas Concentration 

C 
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Average gas concentration inside the pore 


C,* - 




Average Adsorbed Phase concentration 


c; 


C^v 


c. 


Time 


/ =-, r 


L 


•hed 


Distance in radial direction 


^ R-R, 


Introducing the dimensionless variables {Cj,Cp ,Cs*, f and in governing 

equations (3.1), (3.3) and (3.1.3) the following three non-dimensionalised equations 
are obtained: 


+ 


dCr 


+ • 


B dCr. 


D 


1 d 


, +— (^ + ^)(C ‘-C/*) = -— 


I d 


(3.17) 


Here, A 




i?2 i?, 


respectively. 

^ Q\ 1 


B 


\27rV j 


{R,-R,f 

velocity of the gas respectively. 

^ 


Ri and R 2 are the inner and outer radius of the bed 


Q and V are the volumetric flow rate and the superficial 


C= 


(K-R,) 


“'bed 

V V j 


; Lbed is the length of the bed 
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5c /,> rr//") * . r-iSCs 


dC 

here, £ 


^(q;-c/)+f- 


df 


0 


(3.18) 


R 


(k ^ 

(L 

\ 


bed 


yCL j 

V 
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; Rf is the radius of the fiber 


jr;^ f ^ .v max 

a c;,/. 


dC/ 

(T 


Gc;{\~Csy -Hic:) 


N 


(3.19) 


here, G = and H = 

V V 

^.vmax^ ^.smax^ 

Here A, B, C, D, E, F, G and H are intermediate coefficients of the governing 
equations. 

C,' In Eq. (3.18) can also be written as follows. 


C7 \ 


D 


(3.20) 
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The dimensionless initial and boundary conditions can be rewritten as follows, 


at/* = 0, C* 


C' 


a 


c 

'^Gf 


c; = 


(RCy 


l/N 


1 + {KCJ 


IIN 


— c 

r * — 
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c, 


Gf 


(3.21a-c) 


at| = 0, Q* = l, c; = 




,c;=i 


(3.22a-c) 


at^ = l. ^ = 0, ^ = 0,^ = 0 

a^ 


(3.23a-c) 


Numerical solution technique 

For a system of one adsorbing component in an inert gas, we have three 
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partial differential equations (3.1), (3.3), (3.13) containing three dependent variables 
Cg, C p, C.S , which are function of time and radial location. These equations are 
solved simultaneously by finite difference method using a Fortran subroutine D03pcf 
of NAG Fortran library. The Fortran codes are given in the appendix. 



CHAPTER 4 

EXPERIMENTAL STUDIES 

In this section details of the experimental set-up used to study adsorption/desorption 
characteristics of VOCs on ACF are described. The experiments were carried on two 
different types of ACF samples having different BET surface area (one procured 
locally and the other imported) under several operating conditions, such as bed 
temperature, inlet concentration of VOCs, gas flow rate and weight of ACF. Both of 
the above mentioned ACF samples are obtained from viscous rayon precursor. 

4.1 Experimental Set up: 

Figure 4. 1 is the schematic diagram of the experimental set-up used in this study. 
The experimental set up may be assumed to consist of three sections as shown in 
Figure 4.1: (a) gas preparation section, (b) test and (c) analysis sections. In the gas 
preparation section, nitrogen gas was bubbled in a vertical bottle containing the VOC 
liquid at a constant temperature to prepare the gas of desired VOC concentration. The 
bubbler used to bubble the gas was a SS made tube whose bottom end was closed and 
the outer surface were perforated up to a distance of 10 cm from the bottom end, with 
perforated holes of diameter 0.08 mm. The entire perforated section of the bubbler 
was submerged in the VOC liquid. The level of the liquid in the bottle was maintained 
high enough to provide sufficient residence time to the carrier gas for complete 
saturation with VOC. The bottle was kept inside the Freon (R-11) refrigeration unit. 
Therefore, by controlling the temperature the vapor pressure of VOCs and hence, the 
required partial pressure were achieved. Hence, this way the required VOC 
concentration was maintained. The required flow rates of nitrogen over a range of 
0.2- 1.0 slpm were controlled and monitored by mass flow controller (Bronkhorst, 
Netherlands). Prior to bubbling the carrier gas (nitrogen in this case) through the VOC 
liquid, it was passed through two consecutive gas purification sections. The first 
purification section containing silica gel was used to remove moisture and the second 
one containing molecular sieves was used to remove trace amount of hydrocarbons 
which may be present in the carrier gas. The test section consisted of a vertical Teflon 
made reactor (I.D ^ 2.5 cm, O.D == 2.8 cm, L = 10 cm) encapsulated in a SS shell (ID 
^ 4.0 cm, L = 20 cm) with provisions for gas inlet and outlet. Inside the SS shell the 
teflon reactor was vertically and co-axially mounted, with one end closed. The outer 



surface of the reactor was perforated with the holes of diameter 0. 1 mm and center to 
center distance of 0.5 cm. The experiments were carried out on ACF in cloth form 
wrapped over the reactor. Thermocouple used to monitor the bed temperature was 
mounted vertically at the center of the reactor. There was a provision for varying the 
bed temperature from 35 to 100^ C with the aid of electrical heating and PID 
controller (FUJI Electric, Japan). The effluent gas stream from the reactor was passed 
to the analytical section. The analytical section consisted of gas chromatography (GC) 
with FID and data station. A computer was connected to the data station to store the 
chromatograph and the peak areas. 

4.2 Experimental Procedure 

The required amount of ACF sample (~5 g) was weighed by electronic balance 
and then pretreated before can 7 ing out the test runs. The pretreatment was carried out 
at a temperature of ~150 ‘’C and under vacuum for 4 hours. This combination of 
temperature and vacuum pressure were obtained after a series of test runs for the 
optimum performance of ACF. The ACF sample was then wrapped over the Teflon 
reactor. Prior to start of the experiment, the bottle was filled with 140 ml of VOC and 
kept inside the refrigeration unit. For a particular VOC inlet concentration, the 
corresponding saturation temperature of the refrigeration unit and hence the VOC 
bottle temperature was set. Then the refrigeration unit was allowed sufficient time 
(~1.5 hrs.) to attain the required temperature. The reactor was heated to the desired 
bed temperature (-50 ^C) and further kept at that temperature for 1 hr so that the 
system was stabilized and a uniform temperature existed in the bed, The required flow 
rate of the gas was controlled using MFC. The concentration of the inlet gaseous 
mixture in the line bypassing the reactor was measured by GC prior to the start of the 
adsorption process. As the adsorption was started, the concentration of exit gas from 
the reactor was monitored and measured by GC. 

4.3 Product Analysis 

The product gas was analyzed by GC using Porapak and Chromosorb columns 
{length 15 mm, ID = 2.0 mm). The first one was used as a reference column and the 
later one was used to separate nitrogen and toluene gas. Flame ionization detector 
(FID) was used to detect the components at different retention times. The temperature 
of the injector, detector, oven and oven final were maintained constant at 170 °C, 180 



C, 140 C and 150 ^C, respectively. Nitrogen was used as a carrier gas with a flow 
rate of 27 ml/min and a pressure of 3.8 kg/cm^. All the GC data were stored in a 
computer using a data station (Winchrom software) connected between the 
chromatograph and the computer. The areas obtained in a computer were in 
proportion to the concentrations as pre-defined in the software. 

4.4 Experimental variables: 

The experiments were performed in a fixed bed reactor under different operating 
conditions shown in Table 4.1. Most of the possible combinations of bed temperature, 
inlet VOC concentration, gas flow rate, amount of ACF were selected for study. 


Experimental Conditions 

Inlet VOC 
Concentration 
(ppm) 

2000 

4000 

6000 

10000 

Gas flow rate 
(slpm) 

0.2 

0.5 

1.0 

Bed temperature 
(*C) 

35 

40 

50 

_J 

75 

100 

Sorbents amount 

(gm) 

— 

2,2 

4.5 

8,5 

r 

BET areas 

(mVg) 

1000 

(procured locally) 

1700 

(hnported) 


Table 4.1 : Different experimental variables for VOC adsorption 

The experiments were carried out at various VOC inlet concentrations ranging from 
2000 to 10000 ppm, whereas the gas flow rates ranged between 0.2 to 1 .0 slpm. A few 




numbers of experiments were also carried out at six different bed temperatures in the 
range of 35 to 100 °C to study the effect of bpd temperature. Effects of weight of ACF 
on the adsorption of VOCs were also studied for various amount of ACF taken. Here, 
2.2 g of ACF corresponds to 1 wrap, whereas 4.5 and 8.5 corresponds to 2 and 3 
wraps respectively. 

4.5 Regeneration of ACF: 

The schematic diagram for the regeneration (or desorption) of ACF is 
shown in Figure 4.2. The regenerations were carried out to make the ACF available 
for subsequent use. The regenerations of ACF were carried out by electrical (~ DC 20 
V and 3 A) heating under a low nitrogen flow rate of 0.2 slpm. A temperature of ~150 
‘’C and regeneration time of 60-75 min were typically required for the complete 
regeneration of ACF. The temperature was monitored by a thermocouple fixed at the 
surface of the ACF sample. The voltage and current were measured by voltmeter and 
ammeter, respectively. The concentration of the exit gas from the reactor was 
measured with the help of GC using FID. 
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Fig. 4.2: Schematic diagram for the regeneration of ACF 









CHAPTER 5 

RESULTS AND DISCUSSION 

This chapter describes the experimental and model predicted results obtained for the 
removal of VOC by adsorption on ACF under varying operating conditions, such as 
inlet concentration of VOC, bed temperature, flow rate of the gas and weight of ACF. 
The results have been analyzed from breakthrough point of view. The experimental 
results are compared and explained, wherever possible, with the model simulated 
result. 

The term breakthrough refers to the response of an initially clean bed (i.e. free of 
adsorbate) to an influent containing an adsorbate or a mixture of adsorbates and a 
breakthrough curve is the history of the effluent concentration from a bed subject to a 
stepped concentration in input. The shape of a breakthrough curve is determined by 
the nature or type of the adsorption isotherm. In addition, it is influenced by the 
individual transport processes in the bed and within the particles or fibers (macro or 
micropores). In the typical experiments performed to study the effects of several 
operating conditions, the teflon made reactor bed wrapped with fresh ACF was heated 
to a certain bed temperature and fed with a gas containing the calibrated mixture of 
nitrogen and toluene at a constant flow rate. During the initial stage of adsorption 
process, most of the VOC entering the bed is removed by adsorption occurring at the 
surface of the pores. During this stage of adsorption, no or very low concentration of 
VOC is observed at the reactor outlet. However, as the adsorption proceeds with time, 
the adsorption sites of the ACF are occupied by the VOC molecules making them 
unavailable for further adsorption. This phenomenon (saturation of ACF bed) initially 
takes place at the entrance region of the reactor and gradually proceeds towards the 
exit. Thus, in the saturated bed region there is an elutriation and mass transfer is 
confined within a small region (zone). This process continues till the overall 
adsorption process is completed and entire bed is saturated with VOC. Under this 
condition, the VOC concentration at the reactor outlet becomes equal to that at the 
inlet. In the present study, the change in VOC concentration with time at the reactor 
outlet was determined using gas chromatography. Such transient behavior of 
concentration is called breakthrough and the corresponding concentration vs. time 
plot is called breakthrough curve. 
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Various quantities required in the model analysis such as, superficial gas 
velocity, bed porosity, density of the gas, Reynolds number, Schmidt number, 
dispersion coefficient, Peclet number and mass transfer coefficient, were calculated 
under the experimental conditions. These calculations have been presented in the 
appendix. 

5.1 Effect of bed Temperature on adsorbent 

To determine the effects of bed temperature on the VOC removal, the 
experiments were carried out at varying bed temperatures (35, 40, 50, 75 and 100 "C). 
The inlet concentration of VOC (toluene) was kept constant at 10000 ppm (the 
corresponding cooling bath temperature to obtain the desired vapor pressure of VQC 
and hence, the required partial pressure in the gaseous mixture was kept constant at - 
1 .69 *’C). In all the runs, the experiments were carried out on 5 g of the ACF sample 
and at a constant gas flow rate of 0.5 slpm. Figure 5.1 describes the experimentally 
obtained breakthrough characteristics of VOC during dynamic adsorption by ACF 
(Type-1 sample) carried out at five different temperatures. It may be observed from 
Figure 5.1 that the breakthrough times are approximately 7 min for both the 
temperatures: 35 and 40 °C. However, the total adsorption time (time for the exit gas 
concentration to approximately reach the inlet concentration, 1 0000 ppm) is slightly 
greater at 40 ®C (40 min) than at 35 ®C (30 min). However, both the breakthrough and 
adsorption times were found to increase to 10 and 50 min, respectively at the bed 
temperature of 50 '^C. Further increase in bed temperature to 75 ®C resulted in 
significance decrease in both breakthrough time (6 min) and total adsorption time (30 
min). As also observed from the figure, the breakthrough curve further shifted to the 
left as the bed temperature was increased to 100 °C, indicating an early saturation of 
the bed. The similar experiments were carried out on ACF (Type-2 sample), under the 
identical operating conditions to study the effect of bed temperature on the VOC 
adsorption process. The breakthrough curves for the adsorption of VOC on the ACF 
(Type-2 sample) under various temperatures are described in Figure 5.2. It may be 
observed from Figure 5.2 that there was practically no improvement in the 
breakthrough time with an increase in the bed temperature from 35 to 50 °C. In each 
case the breakthrough occurred in approximately 20 min. However the total 
adsorption times were 50, 60 and 70 min for the bed temperatures of 35, 40 and 50 ®C 
respectively. Here also, further increase in bed temperature to 75 °C resulted in 
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decrease in both the breakthrough time (12 min) and the total adsorption time (42 
min). It is also observed from Figure 5.2 that the breakthrough curve shifted sharply 
to the left as the bed temperature was increased to 100 Thus, from these results it 
was concluded that a bed temperature of 50 is most favorable for capturing toluene 
vapor from the effluent and therefore, henceforth all the remaining experiments in this 
study were carried out at 50°C. The breakthrough and total adsorption times for two 
types of ACFs at different bed temperatures are summarized in Table 5.1. The 
difference in the performance of the two types of ACF samples in capturing VOC due 
to the difference in BET surface area and pore size are discussed later in section 5.6. 

The model predictions were done under identical operating conditions to 
explain the breakthrough of VOC in the ACF bed. The VOC adsorption by ACF in a 
fixed bed under isothermal condition is considered to be controlled by three rate 
mechanism as discussed in the theoretical section, 

(a) diffusion of the species from bulk gas to the fiber surface through a gas film 
surrounding the fiber (bulk diffusion). 

(b) diffusion of the gas within the pores (pore diffusion) 

(c) adsorption /desorption of VOCs inside the pores (surface reaction). 

The diffusion of the gas from the bulk to the fiber is dependent on mass transfer 
coefficient, which in turn is determined by the diameter of the fiber and gas flow rate. 
The mass transfer coefficient is shown to be dependent on Sherwood number which in 
turn, is a function of the particles Reynolds number and Schmidt number. Under the 
existing experimental conditions (over a temperature range of 30 to 50 °C) the 
Reynolds number was calculated to be constant at 4.142 x 10'% since in all the cases 
the gas flow rate and the diameter of the fiber were kept constant, whereas Schmidt 
number was found to marginally vary (5%) in the range of 1.38 to 1.45. Hence, in 
each case the effect of bulk diffusion was the same. The diffusion of the gas in the 
pores is primarily dependent on the pore size of the sorbents. With increase in the 
temperature the pore diffusion coefficient is expected to be higher and favorable for 
pore diffusion. However, over a temperature range of 35 to 50 a marginal increase 
in pore diffusion from 1.118x10'^ to 1.145x10'’ m^/s was observed. The main 
reason for longer breakthrough time at higher temperature was attributed to a higher 
rate constant because of its Arrhenius type of temperature dependence. 

If we recall Equations (3.3-3.4) derived in the chapter on theoretical analysis, the 
Sips isotherm equation was obtained by equating the adsorption and desorption rates 
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under equilibrium condition. Therefore, the Sips isotherm parameters b andC^^ have 

been used wherever necessary as adjustable parameters in the present model to 
explain the experimental data. The numerical values of these parameters for the 
propane-activated carbon system have been reported by Do (1998). For the present 
work (toluene-ACF system) these numbers were used as a guideline for explaining the 
breakthrough characteristics of VOC by model prediction. Figures (5.3-5.4) compare 
the model predicted breakthrough curves with the experimental data for the Type-1 
and Type-2 ACF samples, respectively. A reasonable good agreement is observed 
between the two results within the experimental and computational errors. In each 
case the data could be explained by the model with a marginal change in the model 
parameters, b (from 10.0 to 8.5) and (from 0.75 to 0.95) over a temperature range 

of 35 to 50 ”C. The temperature dependency of the Sips isotherm parameters are as 
follows (Do, 1998): 


b = exp 


_Q_ 

R,T 


and C^=C^,oexp 



f, rV 

X 

1- — 




Here, b^ is the adsorption affinity constant at infinite temperature and C is the 


saturation capacity at the reference temperature To. Hence, from the above equations 
it is clear that, with the increase in temperature b will decrease, whereas there will be 
an increase in the values of the parameter as obtained in the present study. The 


adsorption equilibrium constants K, is defined as the ratio of adsorption rate constant 
(ka) to desorption rate constant (ka). The values of ka can be calculated by using the 
formula which is obtained from kinetic theoiy of gases. Once the value of k* is 
obtained, the kd can be calculated by using the formula incorporating ka and which is 
obtained by proper mathematical derivation. The values of K were determined to be 
25.60, 23.94 and 22.82 mVmol for Typel ACF at 35, 40, and 50 ‘’C, respectively and 
21.76, 20.95 and 20.14 mVmol, respectively for Type-2 ACF under the identical 
temperatures. The numerical values of K for two types of samples xmder various 
temperatures are reported in Table 5.2. The exothermic heats of adsorption were 
calculated by using the vant Hoff equation: 




AH 
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Under the existing experimental conditions the heats of adsorptions determined from 
the slope of InK vs 1/T curve were found to be 1.45 and 1.01 kc/mol for Type-1 and 



Type-2 ACF, respectively, as shown in the inset of Figures 5.3 and 5.4. The adjusted 
parameters of Sips isotherm for toluene are reported in Table 5.3 (for Typel ACF) 
and Table 5.4 (for Type 2 ACF). 

5.2 Effect of concentration: 

To determine the effects of VOC concentration at trace levels on the 
breakthrough characteristics, the experiments were carried out for varying VOC inlet 
concentrations: 2000, 4000, 6000, and 10000 ppm. For each run 5 gm of the ACF 
sample was taken. The bed temperature was kept constant at 50 ‘’C and the flow rate 
at 0.5 slpm. Figures (5.5-5.6) describe the experimentally obtained toluene 
breakthrough curves under various inlet concentrations for Type-1 and Type-2 ACF 
samples, respectively. As observed from Figure 5.5, the breakthrough times are less 
than 10 min in each case for Type-1 sample, whereas the total adsorption time 
decreases from 70 to 40 min as the inlet concentration is increased from 2000 to 
10000 ppm. For Type-2 sample, referring Figure 5.6, the breakthrough times are 
observed to be approximately 10 min in each case. However, the total adsorption time 
decreases from 80 to 50 min as the concentration level is increased from 2000 to 
1 0000 ppm. The breakthrough and the total adsorption times are listed in Table 5.5 for 
both types of ACF samples. As also observed from these results, decrease in total 
adsorption time with the increase in concentration is relatively larger at higher 
concentration levels (4000-10000) than lower concentration levels (2000-4000) in 
case of both types of the samples, indicating suitability of ACF in capturing toluene 
vapor at the concentration levels less than 1% (10000 ppm). The increase in the 
breakthrough and total adsorption times with the decrease in inlet concentration levels 
as observed from Figures (5.5-5.6) can be explained in terms of the total amount of 
VOC. With decrease in the inlet concentration under identical flow rates, the total 
amount of VOC (moles) entering the micro-pores of the fiber is less. Therefore, the 
saturation of the ACF bed is delayed and occurs in relatively longer time. 

Figures (5.5-5.6) show that the model predictions agree well with the 
experimental results. For the model predictions of the breakthrough curves under 
identical experimental conditions, variations in the model parameters b and were 
required over the inlet concentration range of 2000 to 10000 ppm. In the case .of 
Type-1 sample b was varied from 23 to 8.5 (kPa)’', while was varied from 0.38 to 

0.95 (mmol/g). Similarly, in the case of type-2 sample b was varied from 45 to 7.5 
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(kPa) ' and C^from 0.49 to 1.70 (mmol/g) for an inlet concentration of 2000 and 

10000 ppm, respectively. In each case, a substantial variation in the calculated values 
of the adsorption equilibrium constant was observed. The value of K in the case of 
Type-1 sample was found to increase from 25.6 to 61.79 m^/mol corresponding to a 
concentration change from 10000 to 2000 ppm. On the other hand, for Type-2 sample 
these values were found to be 20.16 and 120.8 mVmol under the identical 
concentration levels. The slopes of K vs concentration plots as depicted in the 
respective insets of Figures (5.5-5 .6) indicate a nonlinear type of adsorption isotherm 
over the selected concentration range. This is consistent with the Sips isotherm 
assumed in this study for adsorption of toluene on ACF surface. The values of K at 
various concentration levels for both types of ACF samples are reported in Table 5.6. 
As observed from Table 5.6 the adsorption equilibrium constant K in the case of 
Type-2 sample is consistently greater than that in the case of type-1 sample over the 
entire concentration range. Hence, from these results it can be concluded that, the 
adsorption capacity of Type-2 sample is greater than that for Type-1 sample under 
identical concentration levels. The superior performance of the Type-2 sample over 
the Type-1 sample in capturing VOC may be attributed to the difference in BET 
surface area and pore size, which would be discussed later in section 5.6. 

5.3 Effect of flow rate: 

To observe the effect of flow rate on the VOC removal, the experiments were carried 
out at different gas flow rates: 0.25, 0.5 and 1.0 slpm. In each case the bed 
temperature was fixed at the optimum reaction temperature of 50 °C. The weight of 
the sample and the inlet gas concentration were kept constant at 5 gm and 10000 ppm, 
respectively. Figure 5.7 describes the experimentally obtained VOC breakthrough 
curves under the above mentioned flow rates. As observed from Figure 5.7, the 
breakthrough time decreases from 20 to 3 min as the flow rate increases from 0.25 to 
1 slpm. The corresponding times to reach 10000 ppm were 55 and 14 minutes, 
respectively. The breakthrough time and the total adsorption time for all the flow rates 
are listed in Table 5.7. 

Figure 5.7 also compares the experimental data with the model predictions. Good 
agreements are observed between the two results within the experimental and 
computational errors. However, at lower flow rate an overshoot is observed in the 
model predicted breakthrough curve, which is attributed due to the computational 
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error involved in the numerical scheme adopted in this work. A finite difference 
method as per the NAG Foitam subroutine D03pcf was used to solve the partial 
differential equations. In principle, values of the model parameters b and should 

not vary with the variation in the flow rates. Expectedly, as observed fi'om Table 6 the 
values of the adjustable parameters b and were kept constant at 8.5 kPa'^ and 0.95 

mmol/g, respectively for each case. Further experiments were carried out at the 
identical flow rate to determine the repeatability in the data. Figure 5.8 reveals a 
reasonably good repeatability between the data obtained under two sets of 
experiments. 

As also observed from Figure 5.7, both the breakthrough time and the total 
adsorption time decreased with increase in flow rate. This can be explained in terms 
of the particle mass transfer coefficient, k„,. The diffusion of the gas from the bulk to 
the fiber surface is dependent on km, which in turn is a function of Sherwood number, 
Sh, where, _ k^D^ and Dm and Dp are the molecular diffusivity and particle 

* 3/2 — 

A 

diameter, respectively. Sherwood number was calculated fi'om the following 
correlation: Sh = 2.0+1. l(S'c)° (Re)“ ® (Yang, 1 999) 

Here, Re, Reynolds number is dependent on the diameter of the fiber and gas 
superficial velocity. It may be pointed out that in the absence of any correlation 
reported in the open literature for the dispersion coefficient in a radial flow as 
prevalent in the present situation, the above correlation available for determining the 
axial dispersion coefficient in a fixed bed filled with the adsorbent materials was used 
to determine the radial dispersion coefficient over a very small thickness (0.13 cm) of 
the fabric material (ACF). The effects of dispersion coefficients on the breakthrough 
curve will be discussed later in the section on the model parametric study. 

The diaipeter of the ACF sample used was very small (4,6 x 10'!^ m in this 
case). Similarly, the gas velocity in the existing experimental conditions was also 
small in the range between 4.5 x 10*^ and 2.47 x 10'^ m/s, corresponding to flow rates 
of 0.25 slpm and 1.0 slpm, respectively. Hence, due to small fiber size and small gas 
flow rate, the Reynolds numbers were calculated to be very small (almost constant at 
4.14 X 10'^). Therefore, it is evident fi-om the above correlation, the Sherwood 
number became independent of superficial gas velocity {Re) resulting in a constant 
mass transfer coefficient and dependent only on molecular diffusivity and fiber 


diameter. The fiber diameter being very small, the mass transfer coefficient was 
calculated to be very large at 41.24 m/s. The numerical values of various parameters 
used in this model such as superficial velocity, Reynolds number, Sherwood number 
and mass transfer coefficient are reported in Table 5.9. From these results it can be 
concluded that the effect of mass transfer coefficient on the saturation time and hence, 
breakthrough time was negligible under the existing operating conditions. Similarly, 
the pore diffusion coefficient was also calculated to be approximately constant at 1.14 
X 10 m /s over the entire range of flow rate. Therefore, decrease in breakthrough 
time and the total adsorption time with the increase in flow rate is attributed to the 
kinetic rate, which was controlled by the total amount of VOC reaching the surface of 
the fiber, as also found in the case of concentration effects. Table 5.9 also lists the 
numerical values of Schmidt number, dispersion coefficient and effective diffusivity, 
etc. 

5.4 Effect of sorbent amount: 

To study the effect of weight of the adsorbing materials (ACF), the 
experiments were carried out at three different amounts of ACF: 3, 5 and 9 g. These 
amounts of ACF sample corresponds to 1 , 2 and 3 number of wraps, respectively. 
The inlet gas concentration and flow rate were kept constant at 10000 ppm and 0.5 
slpm, respectively. The bed temperature was maintained at 50 °C. 

Figure 5.9 describes the VOC breakthrough data for the various amounts 
ACF sample. As observed from the figure, the breakthrough time increases from 5 to 
16 min with increase in the weight of ACF from 5 to 9 g. The corresponding total 
adsorption times were observed to be 40 to 60 min, respectively. As also observed 
from figure 5.9, the breakthrough came almost instantly for 3 g of the sample, the 
total adsorption time being 30 min. The increase in the breakthrough time with 
increase in the sample amount is obviously due to increase in the total adsorption 
surface area for VOC adsorption. The breakthrough time and the total adsorption time 
for three different amounts of ACF samples are reported in Table 5.8. 

Figure 5.9 shows that the model predictions agree well with the 
experimental results. However, at higher amount of ACF sample an overshoot is 
observed in the model predicted breakthrough curve, as also seen in the case of flow 
rate (Figure 5.7), which is attributed due to the computational error involved in 
numerical system used in the present work. As observed from Table 5.3 the values of 
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the adjustable parameters b and C.^were kept constant at 8.5 kPa ' and 0.95 mmol/g, 
respectively for each case. However, at higher amount of ACF sample (9 g) a lower 
value of Cj^(0.57 mmol/g) was required to explain the experimental data. This may 

again be attributed due to a poor choice of the correlation used in this work to 
characterize the radial dispersion coefficient over a very small thickness (0.13 cm) of 
the bed, as also described in the previous section. 

5.5 Effect of BET surface area: 

The previous breakthrough data obtained for the two types of ACF samples (Type- 
1 and Type-2) under different temperatures and concentrations consistently revealed 
that the breakthrough times and the total adsorption times for Type-2 sample having a 
BET surface area 1700 mVg were greater than that for Type-1 sample having a BET 
area of 1000 m^/g. The superiority of type2 sample over Typel is clearly attributed 
due to a higher specific surface area, responsible for gas adsorption. To elucidate the 
above observation the data are replotted in Figures (5.10-5.11), on a comparable 
format. Under these conditions, the gas flow rate was held constant at 0.5 slpm. In 
either case, the reactor was wrapped with 5 g of the adsorbent. Figures (5.10-5.1 1) 
compares the experimentally obtained breakthrough curves of toluene during 
adsorption by two types of ACF samples (Type-1 and Type-2). As observed from 
Figure 5.10, the breakthrough times and the total adsorption times are greater in case 
of Type-2 sample than those in the case of Type-1 under two different bed 
temperatures (30 and 50 °C). It is clear from Figure 5.10 that at the temperature of 50 
°C the breakthrough time in the case of Type-1 ACF was less than 5 min in 
comparison with 20 min for the Type-2 ACF. Similarly, the times to reach the VOC 
inlet concentration were approximately 40 and 72 min, respectively. As observed 
from Figure 5.1 1, the breakthrough curves follow the same trend under two different 
VOC inlet concentrations (2000 and 10000 ppm), indicating superior performance of 
the ACF with higher BET surface area than the ACF with relatively lower BET 
surface area. 

5.6 Comparative adsorption performance of packed bed and perforated bed 
reactor: 

Figure 5.12 compares the experimental breakthrough curves of toluene obtained from 
packed bed and perforated bed reactors during adsorption by ACF under identical 
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operating conditions. In the former, the ACF samples were packed in the bed after 
pretreatment, while in the later the ACF samples were wrapped over perforated teflon 
made reactor. In each case the experiments were carried out on 5 g of the ACF 
sample and at a constant gas flow rate of 0.5 slpm. The temperature of the bed in both 
cases was maintained constant at 50 °C. At the inlet concentration of 10000 ppm the 
breakthrough time in the case of packed bed was observed to be significantly higher 
(40 min) in comparison with that in the case of perforated bed (15 min), 
corresponding to a total adsorption time of 90 and 50 min, respectively. The same 
trend was also observed for a lower inlet VOC concentration i.e 4000 ppm. Hence, 
from the above observation it can be concluded that the adsorption by ACF in the 
packed form is much more effective than that by the ACF wrapped over the 
perforated teflon reactor under the identical operating conditions. The superiority of 
the ACF sample in packed form relative with that wrapped over the perforated reactor 
may be attributed due to sufficient time of contact between VOC in the laden gas and 
the adsorbent materials in the case of former than later. In the case of packed bed, 

L bed thickness 

residence time is — which is much larger than in the case of 

V V 

*^0 .. 'r 

perforated bed reactor. However, it is due to the ease of regeneration of 
saturated/equilibrated ACF by electrical heating on the ACF wrapped over the reactor 
that the later type of arrangement is generally preferred. 

5.7 Comparative adsorption performance of granular activated carbon (GAC) 
and ACF: 

Figure 5.13 compares the experimental breakthrough curves of toluene during adsorption 
by ACF and GAC under identical operating conditions. In either case, the reactor was 
packed with 10 g of the adsorbents (ACF or GAC). Due to difference in their bulk density 
and therefore, in the length of the bed, the bed porosity was calculated to be 0.68 and 0.37 
for ACF and GAC, respectively. The temperature of the packed bed was held constant at 
45 ° C. The toluene concentration in the inlet gas was 4000 ppm, which corresponded to 
approximately -10 ° C of the bottle temperature. As observed from the figure, the 
breakthrough times in the case of AFC are larger than those in the case of GAC under 
three N 2 flow rates (0.5, 1.0 and 1.5 slpm), indicating superior performance of the former 
than the latter. At higher flow rate (i.e. 1.5 slpm), the breakthrough time in the case of 
GAC was observed to be less than 5 min in comparison with 35 min for ACF. Similarly, 
the times to asymptotically reach the VOC inlet concentration were jqiproximately 5 and 



comparison with GAC is also reported elsewhere (Inagaki, 2000). The superiority of 
ACF over GAC in capturing VOC is attributed to a large number of micropores on their 
surfaces, responsible for the gas adsorption. The model developed in this study was also 
applied to simulate the experimental breakthrough curves. Figure 5.13 compares the 
model results with the adsorption data of toluene for both ACF and GAC. The model 
predictions are observed to be in good agreement with the data. As pointed out in the 
theoretical analysis the values of the adsorption/desorption constants required in the 
model were obtained from a previous work and the remaining parameters were calculated 
based on the existing operating conditions. 

The adsorption experiments were also carried out at relatively higher 
concentration, i.e. 50000 ppm (5%) to determine the relative performances of ACF and 
GAC in adsorbing the organic pollutant. This concentration corresponds to the Jiquid 
toluene temperature of 30 'C. As observed from the breakthrough curves drawn on Figure 
5.14 for varying Nt flowrates (0.2, 0.5, and 1 slpm) the breakthrough times for ACF were 
higher than the respective times for GAC under identical experimental operating 
conditions. The similar breakthrough characteristics were consistently observed under 
several other experimental conditions used in this study. Based on these results, it was 
concluded that ACF exhibited enhanced adsorption capacity for VOC in comparison 
with GAC for the same amount of the materials. 

5.8 Effect of regeneration time on the adsorption performance of ACF: 

In the present study, electrical heating was employed to regenerate the fibers saturated 
with VOC during adsorption. The details of the arrangement were previously described in 
the section on the experiment. A series of the test runs were carried out to determine the 
optimum regeneration conditions, including regeneration time and DC voltage. For each 
case the regenerated fiber cloth was re-subjected to the adsorption test to determine its 
adsorption performance following regeneration. Figure 5.15 describes one of such results 
obtained on the breakthrough of the ACF regenerated for different regeneration times. 
The adsorbent was equilibrated with 8000 ppm of toluene prior to regeneration. The 
voltage was maintained constant at 40 V. The corresponding temperature of the cloth was 
measured to be 1 50 “ C. As observed from the figure, there is insignificant difference 
between the adsorption profiles of the fiber cloth regenerated for 15 and 30 min. The 
adsorption characteristics of the fiber clothes which were regenerated for at least 45 min 
are, however, found to improve, resulting in increase in the breakthrough time. The 
adsorption profiles of ACF regenerated for 60 min are observed to be approximately 
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identical to those of a fresh ACF. Therefore, in this study a regeneration time of 60 
min was chosen to be optimum for the maximum removal of the adsorbed VOC. 
During the regeneration period, the current gradually increased from an initial value 
of 1 A to 3 A, indicating decrease in the electrical resistance of the fiber from 40 Q to 

approximately 13 Q, as the adsorbed VOC was gradually removed from the fabric due 
to heating. 

5.9 Model parametric study: 

Model simulation of a large-scale VOC emission: 

After validation of the model developed in the study, the simulations were carried out for 
an adsorption bed required to control a large scale VOC emission from an industry, A 
reactor of 0.3 m l.D and 1 .0 m height, containing 40 kg of ACF sample wrapped over a 
perforated reactor was assumed. The bed temperature was assumed to be 50 °C and the 
inlet concentration 5000 ppm. Figure 5.16 shows the breakthrough curves for the VOC 
laden flue gas at varying flow rates of 250, 500 and 750 slpm. The corresponding 
breakthrough times are observed to be 300, 200 and 120 min, respectively. Thus with the 
help of above-simulated results, change over time to put a bed of fresh adsorbent 
materials can be estimated. In the present case if the flow rate is 250 slpm a changeover 
time after approximately 5 hr is required. Hence the simulated results are useful in 
designing an industrially sized adsorber. 

Effect of dispersion coefficients: 

Figure 5.17 describes the effect of dispersion coefficient on VOC breakthrough curve 
by model predictions. As discussed in section 5.3 the correlation available for 
determining axial dispersion coefficient in a packed bed was used in the present study 
to characterize the radial dispersion coefficient over a very small thickness (0.13 cm) 
of the fabric bed. This necessitates studying the effects of dispersion coefficient on the 
VOC breakthrough curve under identical operating conditions. The simulations were 
carried out for different dispersion coefficients ranging from 2.92 x 10' to 2.92 x 10’ 
m^/s (difference in five orders of magnitude). As observed from the plots on Figure 
5.17, the breakthrough times and the total adsorption times were approximately same 
in all the cases for which the value of dispersion coefficient was chosen > 2.92 x 10'^ 
m^/s. However, as the dispersion coefficient was decreased to 2.92 x lO"^ m^/s, the 
breakthrough time was observe to increase to 10 min, without any significant change 
in total adsorption time (37 min). With further decrease in the value of dispersion 
coefficient by one order of magnitude (2.92 x 10'^ m^/s) the breakthrough time 
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significantly increased to 20 min and the total adsorption time was found to decrease 
to 28 min. As also observed from Figure 5.17, with further decrease in the dispersion 
coefficient (2.92 x 10 * m^/s) the breakthrough curve became sharp with abrupt rise in 
the concentration level. Beyond this a marginal increase and decrease in the 
breakthrough time and total adsorption, respectively time was observed. 

The above effects (abrupt rise in the concentration level at smaller values 
of the dispersion coefficient) may be attributed due to high inertial effects in the bed. 
As seen from equation 3.17 in the theoretical analysis section the intermediate 
coefficient C involves the dispersion coefficient term. At lower value of dispersion 
coefficient C becomes very small. As a consequence the inertial effects in the bed 
become predominant in comparison with the dispersion effects. Hence, the bed 
appears to be saturated with VOC, with the exit concentration becoming equal to the 
inlet concentration level. 
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Fig-5.1 : Effect of bed temperature on breakthrough of toluene 
(W= 5.0 g, Qxj = 0.5 slpm, C = lOOOOppm) 
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Fig 5.3 : Comparison of model predictions and experimental data for effect of bed temperatur 
determination of heat of adsorption of VOC over ACF-Type 1 
(W= 5.0 g, Qjj = 0.5 slpm, C = 10000 ppm) 
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Fig 5.4 : Comparison of model predictions and experimental data for effect of bed temperature: 
determination of heat of adsorption of VOC over ACF-Type 2 
( W = 5.0 g, = 0.5 slpm, C = 10000 ppm ) 
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Fig. 5.6 : Effect of concentration on breakthrough of toluene over ACF-Type2 
( W = 5.0 g, = 0.5 slpm, T^ed = SO^C ) 
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Fig.5.7 : Effect of gas flowrate on breakthrough of toluene over ACF-Typel 
(W= S.Og, C = 10000 ppm, = 50®C) 
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Fig. 5.8: Effect of gas flowrate on breakthrough of toluene over ACF-Typel 
(W= 5.0g, C = 10000 ppm, 50®C) 
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Fig.5.11 : Effect of BET surface area on breakthrough of toluene over ACF 
under different concentration levels 

(W = 5.0s, Qn, = 0*5slpm, T^ed = SO^^C) 
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Fig. 5.15: Effect of regeneration time on adsorption time by ACF wrapped 

on perforated teflon reactor (Q = 0.2 sipm, W = 6.2 g, = 45 °C) 

(Regeneration conditions: Q = 0.05 sipm, Ti-eri = 150 °C, Volts = 40 v) 





Fig. 5.16: Model predictions for VOC removal on an industrial scale 
(L = 1.0 m, I.D = 0 J m, W = 40 kg, T = 50 ®C, 

Cjniet ~ ppm, thickness = 0.254, m) 
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Fig. 5.17: Effect of dispersion coefficient on VOC breakthrough 
( Rp = 3 X lO'^ m, = 41.247 m/s, = 9.47 x 10"^ m^s) 

(T = SO^C, = 10000 ppm, W = 5 g) 




Sample 

S. No. 

Bed 

temperature 

("C) 

Breakthrough 
time (min) 

Total 

adsorption 
time (min) 

ACF-Type 1 

(Procured locally) 

1 

35 

7 

30 

2 

40 

7 

40 

3 

50 

10 

50 

4 

75 

7 

25 

5 

100 

3 

21 

ACF-Type 2 

(Imported) 

6 

35 

20 

50 

7 

40 

20 

60 

8 

50 

20 

j 

70 

9 

75 

12 

42 

10 

100 

i 

8 

40 


Table 5.1: Breakthrough time and total adsorption time of two different type of 
ACF samples at different bed temperatures. 


Temperature(®C) 

K (m^/mol) 

ACF (Type-1 sample) 

(Procured locally) 

ACF (Type-2 sample) 

(Imported) 

35 

25.60 

21.76 

40 

23.94 

20.95 

50 

22.816 

20.14 


Table 5.2: Adsorption equilibrium constant at various bed temperatures 





Fig. 5.3 : Effect of bed temperature on Breakthrough of Toluene over ACF 
( W = 5 g, Qn 2 = 0.5 slpm, C = 10000 ppm) 

T emperature(^C) 

b (1/kPa) 

( mmol/gm) 

35 

10.0 

0.75 

40 

9.2 

0.85 

50 

8.5 

0.95 

Fig. 5.5 : Effect of concentration on breakthrough of toluene over ACF 
( W = 5 g, Qn 2 = 0.5 slpm = 0.5 slpm, Tbed = 50"C) 

Concentration (ppm) 

b (1/kPa) 

( mmol/gm) 

2000 

23.0 

0.38 

4000 

13.0 

0.64 

6000 

12.0 

0.80 

10000 

8.5 

0.95 

Fig. 5.7 : Effect of gas flowrate on breakthrough of toluene over ACF 
(W = 5 g, C = 10000 ppm, Tb«d = 50“C) 

Flowrate (slpm) 

b (1/kPa) 

( mmol/gm) 

0.25 

8.5 

0.95 

0.50 

8.5 

0.95 

1.00 

8.5 

0.95 

Fig. 5.8: Effect of weight on Breakthrough of toluene over ACF 
(Q = 0.5slpm, C = 10000 ppm, Tbed = 50"C) 

Weight (g.) 

b (1/kFa) 

( mmol/gm) 

3 

8.5 

0.95 

5 

8.5 

0.95 

9 

8.5 

0.57 


Table 5.3: Sips isotherm parameters for various experimental conditions 
(Sample # ACF-Typel, N = 1.0) 



Fig. 5.4 : Effect of bed temperature on Breakthrough of Toluene over ACF 
( W = 5 g, Qn 2 = 0.5 slpm, C = 10000 ppm) 

Temperature(^C) 

b (1/kPa) 

( mmol/gm) 

35 

8.5 

1.33 

42 

8.0 

1.45 

50 

7.5 

1.70 

Fig. 5.6 : Effect of concentration on breakthrough of toluene over ACF 
( W = 5 g, Qn 2 = 0.5 slpm = 0.5 slpm, Tbcd = 50®C) 

Concentration (ppm) 

b (1/kPa) 

c.. 

( mmol/gm) 

2000 

45.0 

0.49 

4000 

25.0 

0.90 

6000 

13.9 

1.14 

lOOOO 

7.5 

1.70 


Table 5.4: Sips isotherm parameters for varioas eiperimental Conditions 
(Sample # ACF-Type2, N = 1.0) 






Sample 

S. No. 

VOC inlet 
concentration 

Breakthrough 
time (min) 

Total 
adsorption 
time (min) 

ACF-Type 1 

(Procured locally) 

1 


5 

70 

2 


6 

65 

3 

6000 

6 

52 

4 

10000 

7 

40 

ACF-Type 2 

(Imported) 

5 

2000 

12 

80 

6 

4000 

12 

62 

7 

6000 

10 

60 

8 

10000 

9 

50 


Table 5.5; Breakthrough time and total adsorption time of two different type of 
ACF samples at different VOC inlet concentration levels. 


Concentration (ppm) 

K(m^/mol) 

ACF (Type-1 sample) 

(Procured locally) 

ACF (Type-2 sample) 

(Imported) 

2000 

61.79 

120.84 

4000 

34.91 

67.14 

6000 

32.22 

37.32 

10000 

25.60 

20.14 


Table 5.6: Adsorption equilibrium constant at various VOC inlet concentration 


Levels 








S.No 

Flow rate 
(slpm) 

Breakthrough 

Time 

(min) 

Total 

Adsorption 

Time 

(min) 

1 

0.2 

20 

55 

2 

0.5 

5 

45 

3 

1.0 

3 

14 


Table 5.7: Breakthrough time and total adsorption time at various gas flow rates 


S.No 

Adsorbent 

amount 

(gm) 

Breakthrough 

Time 

(min) 

Total 

Adsorption 

Time 

(min) 

1 

1 

3 

~0 

30 

2 

5 

5 

40 

3 

9 

16 

60 


Table 5.8: Breakthrough time and total adsorption time at different amounts of 


adsorbent 








CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 

The various conclusions that can be drawn from this study are summarized as follows: 

1. The experimental results revealed that activated carbon fiber is a potential 
adsorbent for capturing VOCs at low concentration levels, due to its high BET surface 
area and micropore volume. 

2. A bed temperature of 50 °C was found to be most favorable for capturing toluene 
vapor from a nitrogen stream in terms of longer breakthrough and adsorption times. 
The dynamic adsorption experiments carried out below and above this temperature 
resulted in smaller breakthrough and adsorption times. 

3. Under the existing experimental conditions (T~ 35-50 “C) the heats of adsorptions 
were found to be 2.06 and 1.02 kc/mol for Type-1 and Type-2 ACF samples, 
respectively. 

4. The breakthrough time during adsorption was observed to considerably decrease 
with increase in the concentration levels suggesting utility of the adsorption method 
for controlling VOC from gaseous effluent at ppm and sub-ppm levels. For example, 
the total adsorption time was found to increase from 40 to 70 min as the concentration 
level wss decreased from 10000 to 2000 ppm. 

5. The commercially obtained ACF was shown to exhibit greater adsorption for VOC 
than granular activated carbon under identical operating conditions. This was 
attributed due to higher BET surface area of ACF (-1200 m^/g) in comparison with 
that for GAC (-600 m^/g). 

6. Adsorption by ACF in the packed form was found to be more effective than that by 
the ACF wrapped over the perforated teflon reactor under the identical operating 
conditions (equal amount of adsorbent weight, temperature, concentration and flow 
rate). However, the ease of regeneration of saturated/equilibrated ACF by electrical 
heating in the case of ACF wrapped over the reactor was the main reason for using 
such type of arrangement over packed bed reactor. 

7. Under existing operating conditions the ACF sample was regenerated 
approximately 25-30 times between two subsequent adsorption experiments, without 
exhibiting any degradation in its adsorption performance (breakthrough time). 

8. A mathematical modeling for the prediction of breakthrough curve was developed 
based on three governing equations: (1) the species balance of the component in the 
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perforated bed, (2) species balance inside the pores of the fiber, and (3) 
adsorption/desorption rates within the pores. The Sips isotherm was found appropriate 
in explaining adsorption/desorption behavior of toluene on ACF. The experimental 
and model predicted results were found to be in good agreement within the 
experimental and computational error. 

9. The fiber diameter being very small, the mass transfer coefficient in the bed was 
calculated to be very large at 41.24 m/s, even at small particle Reynolds number of 
4.14.x 10'^. This suggests that the resistance to the gas film mass transfer coefficient 
was negligible and the adsorption became independent of the film resistance. 

Scope of future study: 

The following recommendations are made to carry out further studies: 

1. Studies on the removal of some other toxic VOCs (such as benzene, xylene, 
dimethylchloride, and trichloroethylene) and other gaseous pollutants like NOx (NO 
& NO 2 ), SO 2 , CO may be carried out on ACF or with other types of adsorbent 
materials like carbon molecular sieves (CMS), activated charcoal and zeolite. 

2. Separate research studies may be carried out on the physical and chemical 
treatments required in carbonization and activation. of raw fibers to enhance the BET 
surface area, pore size and to modify the surface functional groups which are 
responsible for the adsorption of various gaseous pollutants. 

3. On theoretical sides, efforts can be made to find out a suitable correlation for the 
radial dispersion coefficient to achieve more accurate model predictive results. 
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APPENDIX 


The following are the sample calculations for the parameters used in the model. The 
calculations have been done corresponding to the experimental conditions: 


Superficial gas velocity: 
gas flow rate = 0.5 slpm 
Superficial gas velocity = 

where flow rate is in slpm 
Superficial velocity 


(Flowrate)(R)(Temp) 
{22A)(60){pressure){n:DLff^.j ) 

(0.5)(8.314)(323) 


(Al) 


(22.4)(60)(1 0 1 325)(?z- x 0.0254 x 0. 1) 


= 1.23x1 O'-’ m/s 


Residence time: 


Residence time of gas in the reactor = 


Bed thickness 
vel 


Residence time = 


1.27x10"^ 
1.23x1 O’-’ 


1.023 s 


{A2) 


Concentration of VOC in feed stream: 


^ . (pressure)(mol fraction) 

Concentration =— 

{R)(Temperature) 


(101325)(0.01) 

(8.314)(323) 


0.2263 mol/m^ 


(A3) 


Bed porosity 

Amount of sorbent = 5 gm 

Solid density of the sorbent = 1 200 kg/ 

Actual volume of ACF sample = = 4.166x10”^ m^ 

1200 

Bulk volume of the adsorbent bed = ;r[(0.014)^ -(0.0127)‘Jx0.1 = 1.09x10”^ m"* 


Bed porosity e = 


(1.09xl0~^-4.166xl0~^) 

1.09x10"' 


0.647 


(A4) 



Density of gas 


( pressure)(mol. wt) 
{temperature)(R) 


(101325)(92jc10-^) 

(8.314)(323) 


= 1.056 kgW 


(A5) 


Reynolds number based on fiber diameter 

(gas density){ vel)( fiber diameter) 

Kc — 

(gas viscocity) 

(1.056)(1.23xl0'^)(4.6xlQ~^) 

(1.45x10'') 

Schmidt number 

(gas viscosity) 


(A6) 


= 4.14x10 


-5 


& = 


(gas density)(binary molecular diffusivity) 
1.45x10"' 


(A7) 


= 1.45 


(1.056)(9 .47x10'*) 

Dispersion co-efficient 

The radial dispersion coefficient has been approximated by the following correlation 
for axial dispersion coefficient in a fixed bed (Yang, 1997): 


D. 


=^(20 + 0.550 Re) 


(A8) 


9.47x10' 

0.647 


9.47x10' 

0.647 


■ [20 + 0.5 (1.45)(4.14xl0'')] 


-[20 + 0.5 (1.45)(4.14xl0'')] 


= 2.92x10 m/s 


(A9) 


Mass transfer coefficient: 

The mass transfer co-efficient is calculated using the following correlations: 
5/1 = 2.0 + 1. 1(50)°" (Re)®* 



_ (2 + \ A(Sc)'^ {Rey^){moLdiffiismty) 
{fiber diameter) 

= 41.24 m/s 


Calculation of diffusivity: 


The effective diffusivity inside the pores is given by a combination of both Knudsen 
and molecular diffusion. 

The Knudsen diffusivity is given by 


2xr ZRT - T - 

Dk = ^(— =91xR(—Y 

^ 3 nM 


(A 10) 


where 

Rpore= pore radius in m 

M = Molecular weight in gm/gmmol 

T = Temperature in K 

Dk = Knudsen diffusion in m^/s 

The Molecular diffusivity for a binary gas mixture is given by the Wilke-Lee equation 
(Treybal, 1925); 


A, 


1 O'" (1 .084 - 0.24(1 / + 1 / MJ-' )r^ (1 / M^+l/MJ 


\0.5 




kf 

s 


\ 


(All) 


AH ) 


Here, 

Dm= Molecular diffusivity, m^/s 
T = Temperature in K 

Ma, Mb = Molecular weight of the two interacting species (A: Toluene, B: N 2 ) 
Rab = molecular separation at collision, nm = (rA+rB)/2 
e^„= energy of molecular attraction = {s^sff 
k = Boltzmann’s constant 


xO.5 


/ 


^ = collision function, a function of 


\^AH ) 


'AB 


rA = 0.602 m, ^ = 464.398 if 
k 

re = 0.380 m, ^ = 71.4A: 
k 



^ = (464.398x71.4)“' = 182.09 AT 
k 


A 


^ kT^ 


\^AHj 


r-i 


= 0.56 

0.3798 + 0.602 


TolinfiK’— N > 


= 0.4909 nm 


(1/M, + 1/M«)“' =(1/92 + 1/28)“' =0.2158 
Substituting all the values in A(1 1) 


D,„ = 


_ 1 0"' (1 .084 - 0.24(0.2 1 58)° ')(323)''(0.2 158) 


\ 0.5 


= 9.472x10'“ m7s 


(101325)(0.4909)'(0.56) 

The combined diffusivity contributed by both molecular and Knudsen diffusion is 
given by 


_L__L _1_ 

D~ D.., 


D = Diffusion co-efficient 

The effective diffusivity inside the pores is found by using the formula 
D=^ 


(A12) 


(A13) 


where 

a= Intrafiber void fraction 
D = diffusivity co-efficient 
r = tortuosity factor 
For micropore radius of 3.0 x 10'^ m 


The Knudsen diffusivity = (97)( 3.0 x 10'^ ) 


323 V 

. 92 J 


= 5.45 X 10'^ m7s 

The Diffusion co-efficient is given by 

1 11 
D ~ 5.45x10'’ ^9.472x10'“ 



D = 5.15x lO'^Ti'/s 


The effective diffiisivity is calculated as follows 


where r = 5 


_ (5.15xl0'’)x0.7 ,,, 2, 

n =2^ L = 1.14x10 ’ mis 


Adsorption and desorption rate constants: 

From kinetic theory of gases adsoqjtion rate constant, ka is defined as: 

, fifrxioooY' , . , 

fCa - ^ » where s is the sticking coefficient 

V. 2?iM ) 

Here, s = 1.0 

, (^8.314x323x1000^ , 

k=\ =68.15 m/s 

2x;rx92 J 

Sips isotherm is given as follows (Do, 1998), 

1 

(hpY 

^ ^ r ’ P ^ / g 

\ + (bp)^ 

Here Cs can be written as follows, 

\ nJli / ^ ijtLL 

^ ^ l+(^>p)" 

Substituting p = CpRT x lO’^ and comparing with equation (3.4) in theoretical 
analysis section we have , 

^ = bRTx lO ' 

K 

Here Csmax in equation (3.4) is given as follows. 


(A14) 


C =-^xl0" 
BET 


Desorption rate constsn.t,k^ = 


k 68.15 

a — 

bRTxW 12x8.314x323x10" 


2.11mol/s-m" 



25000 


Input files for the code: 

Oper: 

0. 5 0.01 1.0132x10^ 50.0 92.0 1.0 1x10-^ 

1 . Flow rate of Gas (slpm) 

2. Inlet Mol fraction of VOC 

3. Pressure (N/m^) 

4. Initial Temperature (°C) 

5. Molecular Weight, 

6. Stick (sticking probability, to calculate rate constant for adsorption) 

7. Initial mole fraction ofvoc 

8. Final time step, s 

Bed: 

0. 1. 5.0 1200 0.55 0.5 700 

1. Length of bed, m 

2. Weight of the ACF, gm. 

3. Solid density of the ACF, kg/m'' 

4. Outer radius of the bed, inch, 

5. Inner radius of the bed, inch. 

6. Bulk density of ACF, kg/m^ 

Inert 

1.45x10-^ 28.0 

1 . Viscosity of gas, pa-s 

2. Molecular weight of gas 

Sorbent: 

3.0x10'® 5.0 4.6 xlO'"^ 1000.0 

1 . Radius of Pore, m 

2. Tortuosity 

3. Diameter of fiber, m 

4. BET surface area, m^/gm. 

Vocc3: 

8.5 0.95 1.0 9.472x10-* 


toluene 

1. bp 

2. Csv (m.moFgm) 

3. fs 

4. Molecular dififusivity, mVs 



C Last change: I 28 Jan 2003 3:51 pm 

* D03PCF Example Program Text 

* Mark 15 Release. NAG Copyright 1991. 

C Last change: 1 26 May 100 9:31 pm 


* .. Parameters .. 

INTEGER NOUT 
PARAMETER (NOUT=6) 

INTEGER NPDE, NPTS, INTPTS, ITYPE, NEQN, NIW, NWK, NW 
PARAMETER (NPDE=3,NPTS=60000,1NTPTS=6,1TYPE=1 .NEQN=NPDE*NPTS, 
+ N1W=NEQN+24.NWK=(10+6*NPDE)*NEQN, 

+ NW=NWK+(21 +3*NPDE)*NPDE+7‘NPTS+54) 

* .. Scalars in Common .. 

* .. Local Scalars .. 

DOUBLE PRECISION ACC, HX, PIBY2. TOUT. TS, ALPHA 
INTEGER I. IFAIL, IND, IT, ITASK, ITRACE, M 

* .. Local Arrays .. 

DOUBLE PRECISION U(NPDE,NPTS), UOUT(NPDE, INTPTS, ITYPE), W(NW), 

+ X(NPTS), XOUT(INTPTS) 

DOUBLE PRECISION Cgf,cgi,Csmax,adcet,Fmass,outrad,innraci,ftime, 

+ Volflow,Eps,Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid, 

+ Deff,Vmf,Diff,bp,csv,temp,dia,Dz,yi,wfib,denf,res, 

+ Gden,dens,Cpg,Cps,Ya,pi,aprime,Vel,esa,Vo,bden, 

+ Vmw,stick.Rcons, AS,BS,CS,DS,ES,GS,FS,decet,aec,HS 


INTEGER IW(NIW) 

.. External Functions .. 

DOUBLE PRECISION X01AAF 
EXTERNAL X01AAF 

.. External Subroutines .. 

EXTERNAL BNDARY, D03PCF, PDEDEF, UINITID03PZF 

.. Intrinsic Functions .. 

INTRINSIC SIN 
.. Common blocks.. 

COMMON A/BLE/ALPHA 

COMMON/prop/Cgf,cgi,Csmax,Fmass,outrad,innrad,Volflow,Eps, 

1 Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid,aprime, 

2 Deff,Vmf,Diff,bp.csv,temp,dia,Dz,yi,Vel,esa 
COMMON/prop1/Gden,dens,Cpg,Cps 
COMMON/prop2A/mw,stick,pi,Rcons 
COMMON/prop3/decet,aec,adcet 
COMMON/Prop4/AS,BS.CS.DS,ES,GS,FS,HS 

OPEN(UNIT=1,FILE="c:\deba\oper.bct",STATUS="UNKNOWN") 

OPEN(UNlT=2.FlLE="c:\deba\bed.txf,STATUS="UNKNOWN") 

OPEN(UNIT=3,FILE="c:\deba\iner.bct",STATUS="Unknown’') 

OPEN(UNIT=4,FILE=''c:\deba\sorbent.txt".STATUS="UNKNOWN") 

OPEN(UNIT=5,FILE="c:\deba\vocc3.b(t",STATUS="UNKNOWN") 

OPEN(UNIT=8,FILE=’'c:\deba\adsorb1.bd".STATUS="UNKNOWN") 

OPEN(UNIT=9,FILE="c:\deba\adsorb2.txt",STATUS="UNKNOWN") 

OPEN(UNIT=10,FILE="c:\deba\adsorb3.txt",STATUS= 

1 "UNKNOWN") 



OPEN(UNIT= 1 1 .FILE=''c:\deba\adsorb4.txt",STATUS= 

1 "UNKNOWN") 

.. Executable S 

READ(1 ,*) volflow.VMF, Pres, temp, Vmw.stick.yi 
WRITE(1 1 ,*)"volflow,VI\/IF,Pres,temp,Vmw, stick, yi" 

WRITE(1 1 ,*)volflow,VMF,Pres,temp,Vmw,stick,yi 
READ(2,*) Lbed,wfib,denf, outrad,innrad,bden 
WRITE(1 1 ,*)"Lbed,wfib,denf,outrad,innrad,bden " 

WRITE(1 1 ,*) Lbed,wfib,denf,outrad,innrad,bden 
READ(3,*) Gvis,Gmw 
WRITE(11,*)"Gvis,Gmw" 

WRITE(11,*) Gvis,Gmw 
READ(4,*)Rpore,tor,diaf,surarea 
WRITE(1 1 ,*)"Rpore,tor,diaf,surarea" 

WRITE(1 1 ,*)Rpore,tor,diaf,surarea 
READ(5,*) bp,Csv,fa,Diff 
WRITE(11,*) "bp,Csv,fa,Diff" 

WRITE(1 1 ,*)bp,Csv,fa,Diff 
temp=temp+273.0 ! temperature,k 
rcons=8.314 !gas conbstant,(N/m2)(m3)/mol.k 
PI=4.0D0*ATAN(1 .OeO) 
dia=2.0*innrad ! inch 
dia=dia*2.54e-2 !dia of bed, inch to meter 
outrad=outrad*2.54e-2 !dia of bed, inch to meter 
innrad=innrad*2.54e-2 !dia of bed, inch to meter 

Vo=PI*(outrad**2-innrad**2)*Lbed ! Bulk volume of the adsorbent bed, m3 
eps=(Vo-wfib/(denf*1000.0))A/o ! Bed porosity 

pvoid=surarea*rpore*bden*1 000.0/2.0 ! intrafiber porosity, dimensionless 
aprime=surarea*denf*1 000.0 ! Total adsorption surface area per unit volume of fiber,rr 
Cgi=pres*yi/(Rcons*temp) ! Initial cone. moles/m3 

Vel=Volflow*Rcons*Temp/(22.4*60.0*Pres*PrDia*Lbed) Isuperficial gas Velocity,m/s 
Gden=Pres*Gmw*0.001/(Rcons*Temp) IGas density,kg/m3 
Re=Gden*Verdiaf/Gvis ! Reynolds Number 
sc=gvis/(Gden*Diff) Ischmidt number 

Fmass=Diff/diar(2.0+1 .1*(Sc)**0.333*(Re)**0.6) IMass Transfer coefficient,m/s 
pin=pres*VMF Ipartial pressure, pa 
Cgf=Pin/(Rcons*Temp) linletconc. moles/m3 
Dk=97.0*Rpore*(Temp/Vmw)**0.5 IKnudsen diffusivity,m2/s 
Dz=Diff/eps*(20.0+0.5*Sc*Re) laxial dispersion coefficient, mils 

Deff=Dk*(Pvoid/tor) !m2/s, effective diffusivity inside the pores 
adcet=stick*(rcons*1.0e3*Temp/(2.0*PrVmw))**0.5 ladsorption coefficient, m/s 
decet=adcet/(bp*Rcons*Temp*1.0e-3) Idesorption coefficient, moles/s.m2 
aec=adcet/decet ladsorption equilibrium constant 

esa=4.0/diaf ! external surface area per unit volume of fiber,1/m 
Csmax=Csv*1 .0/surarea*1 .Oe-3 Imoles adsorbed/BET area (sqm),Sips equation para 
AS=innrad/(outrad-innrad) Icoefficient of the governing equation 
Ya=Volflow*1 .Oe-3/(2.0*pi*Vei*60.0)!intermediate coefficient of the governing equation 
BS= Ya*(1 .0/(outrad-innrad))**2 Icoefficient of the governing equation 
CS=(Dz/(outrad-innrad)**2)*(Lbed/Vel) Icoefficient of the governing equation 
DS=((1.0-Eps)/Eps)*Fmass*(Lbed/Vel)*esa Icoefficient of the governing equation 
ES= 6 . 0 /(diarPvoid)*Fmass*(Lbed/Vel)!coefficient of the governing equation 
FS=(aprime/Pvoid)*(Csmax/Cgf) Icoefficient of the governing equation 
GS=(adcet*Lbed*Cgf)/(Vel*Csmax) 

HS= (decet*Lbed)/(Vel*Csmax) 



res=LbedA/el ! sec 


write(1 1 ,*) "temp=",temp !temperature,k 
write(1 1 ,*) "clia=",dia !dia of the bed.m 
write(11,*) "PI=",Pl 

write(1 1 ,*) "Cgi=",Cgi iinitial concentration of voc,nnol/m3 

write(1 1 ,*) 'Vel=",vel Isuperficial gas Velocity, m/s 

write(1 1 ,*) ''Gden=",Gden !gas density, kg/m3 

write(1 1 ,*) " Re=",Re IReynolds Number 

write(1 1 ,*) "sc=",sc Ischmidt number 

write(1 1 ,*) "Fmass=",Fmass Imass transfer coefficient, m/s 

write(1 1 ,*) " pin=", pin ! partial pressure, pa 

write(1 1 ,*) "Cgf=",Cgf Ifeed concentration of voc,mol/m3 

write(1 1 ,*) "Dk=",Dk IKnudsen diffusivity,m2/s 

write(11,*) "Deff=",Deff !m2/s, effective diffusivity inside the pores 

write(1 1 ,*) "Dz=",Dz laxial dispersion coefficient 

write(1 1 ,*) "adcet=",adcet ladsorption coefficient,m/s 

write(1 1 ,*) "decet=",decet Idesorption coefficient, moles/s. m2 

write(1 1 ,*) "aec=",aec ladsorption equilibrium constant 

write(1 1 ,*) ''esa=",esa ! external surface area per unit volume of fiber,1/m 

write(1 1,*) "Csmax=",Csmax 

write(1 1 ,*) "AS=",AS 

write(1 1 ,*) "Ya=",Ya . 

write(11,‘) "BS=",BS 

write(11,*) "CS=",CS 

write(11,*) "DS=",DS 

write(1 1 ,*) "ES=",ES 

write(11,*) "FS=",FS 

write(11,*) "GS=",GS 

write(11,*) "HS=",HS 

write(11,*) "eps=",eps IBed porosity 

WRITE(11 ,*)"aprime=",aprime lAdsorption surface area per unit volume of fiber 
WRITE(1 1 ,*)"aipha=",pvoid ! Intrafiber void fraction 
write(11,*) "res=",res 

* .. Data statements .. 

DATA XOUT(1)/0.0e+0/, XOUT(2)/0.40e+0/, 

+ XOUT(3)/0.6e+0/, XOUT(4)/0.8e+0/, 

+ XOUT(5)/0.9e+0/, XOUT(6)/1 .Oe+0/ 


ACC = 1 .Oe-4 
M = 1 

ITRACE =0 
ALPHA = 1.0e0 
, IND = 0 
ITASK = 1 

* 

* Set spatial mesh points 

* 

PIBY2 = 0.5e0*X01 AAF(PI) 
HX = PIBY2/(NPTS-1) 
HX=1.0/(npts-1) 

X(1) = O.OeO 
X(NPTS) = 1 .OeO 
DO I = 2, NPTS - 1 
X(l)= S1N(HX*(I-1)) 



X(l) = HX*(I-1) 
enddo 

★ 

* Set initial conditions 
*• 

TS = O.OeO 
TOUT =0.00001 

WRITE (1 1 ,99999) ACC, ALPHA 
WRITE (1 1 ,99998)(XOUT(l),l=1 ,6) 

★ 

* Set the initial values 

CALL UINIT(U,X,NPTS) 

DO IT = 1,500 

! IF(tout*res.ge.ftime)GO to 1 11 

FAIL = -1 
TOUT=TOUT+0.2 

CALL D03PCF(NPDE,M,TS,TOUT,PDEDEF,BNDARY,U,NPTS,X,ACC,W,NW,IW, 
+ NIW,ITASK,ITRACE,IND,IFAIL) 

* Interpolate at required spatial points 

* CALL D03PZF(NPDE,M,U,NPTS,X,XOUT,INTPTS,ITYPE,UOUT,IFAIL) 
WRITE(10,99996)TOUT*res/60.0,U(1,NPTS/4),U(1 ,NPTS/2), U(1 ,NPTS) 
enddo 

* 

* Print integration statistics 

WRITE (11,99997) IW{1), IW(2), IW(3), IW(5) 

STOP 

99999 FORMAT {jP Accuracy requirement = ',e12.5,/' Parameter ALPHA =', 

+ ’ ',e12.3./) 

99998 FORMAT (' T / X ',6F8.4,/) 

99997 FORMAT (' Number of integration steps in time 

+ \A/ Number of residual evaluations of resulting ODE sys', 

+ tern', 14,/' Number of Jacobian evaluations 
+ ' ',14,/' Number of iterations of nonlinear solve', 

+ 'r ',14,/) 

99996 FORMAT (1x,F14.6,3(1x,F14.6)) 

END 

SUBROUTINE UINIT(U,X,NPTS) 

* Routine for PDE initial conditon 

* .. Scalar Arguments .. 

INTEGER NPTS 

* .. Array Arguments .. 

DOUBLE PRECISION U(3.NPTS), X(NPTS) 

DOUBLE PRECISION Cgf,cgi,Csmax,adcet,Fmass,outrad,innrad, 

+ Volfiow,Eps,Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid, 

+ Deff,Vmf,Diff,bp,csv.temp,dia,Dz,yi, 

+ Gden,dens,Cpg,Cps,Ya,pl,aprime,Vel,esa, 

+ Vmw,stick,Rcons,AS,BS,CS,DS,ES,GS,FS,decet,aec,HS 


.. Scalars in Common .. 



.. Local Scalars .. 

INTEGER I 
.. Scalars in Common .. 

DOUBLE PRECISION ALPHA 

.. Common blocks .. 

COMMON A/BLE/ALPHA 

COMMON/prop/Cgf,cgi,Csmax,Fmass,outrad,innrad,Volflow,Eps, 

1 Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid,aprime, 

2 Deff,Vmf,Diff,bp,csv,temp.dia,Dz,yi,Vel,esa 
COMMON/prop1/Gden,dens,Cpg,Cps 
COMMON/prop2A/mw, stick, pi , Rcons 
COMMON/prop3/decet,aec,adcet 
COMMON/Prop4/AS,BS.CS.DS.ES.GS,FS,HS 


.. Executable Statements .. 

DO 20I=1,NPTS 
U{1.l)=Cgi/Cgf 
U(2,l)=Cgi/Cgf 

U(3.l)=(aec*Cgi)**(1 .0/fa)/((aec*Cgi)**(1 .0/fa)+1 .0) 

20 CONTINUE 

RETURN 

END 

SUBROUTINE PDEDEF{NPDE,T,X,U,DUDX,P,Q,R,IRES) 

.. Scalar Arguments .. 

DOUBLE PRECISION T, X 
INTEGER IRES, NPDE 
.. Array Arguments .. 

DOUBLE PRECISION DUDX(NPDE), P(NPDE,NPDE), Q(NPDE), R(NPDE), 
+ U(NPDE) 

DOUBLE PRECISION Cgf,cgi,Csmax,adcet,Fmass,outrad,innrad, 

+ Volflow,Eps,Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid, 

+ Deff,Vmf,Diff,bp,csv,temp,dia,Dz,yi, 

+ Gden,dens,Cpg,Cps,Ya,pi,aprime,Vel,esa, 

+ Vmw,stick, Rcons, AS,BS,CS,DS,ES,GS,FS,decet,aec,HS 

.. Scalars in Common .. 

DOUBLE PRECISION ALPHA 

.. Common blocks .. 

COMMON A/BLE/ALPHA 

COMMON/prop/Cgf,cgi,Csmax,Fmass,outrad,innrad,Volflow,Eps, 

1 Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid,aprime, 

2 Deff,Vmf,Diff,bp,csv,temp.dia,Dz,yi,Vel,esa 
COMMON/prop1/Gden,dens.Cpg,Cps 
COMMON/prop2/Vmw,stick,pi.Rcons 
COMMON/prop3/decet,aec,adcet 
COMMON/Prop4/AS,BS,CS,DS,ES,GS,FS,HS 



.. Executable Statements .. 


Cfs= (U(1)*Fmass+1 0.0/diaf‘Deff*U(2))/(Fmass+10.0/cliaf*Deff) 
Q(1)=BS/(CS*X)*DUDX(1)+ DS/(CS*X)*(X+AS)*(U(1)-Cfs) 

Q(2)= -ES*(U(1)-Cfs) 

Q(3)= -GS*U(2)*(1 -U(3))**fa+HS*U(3)**fa 
R(1) =(1+AS/X)*DUDX(1) 

R(2) =0.0e+0 
R(3) =0.0e+0 
P(1,1) =(X+AS)/(CS*X) 

P(1,2) = 0.0e+0 
P(1,3) = O.Oe+0 
P(2,1) = O.Oe+0 
P(2,2) =1 .OeO 
P(2,3)=FS 
P(3,1) = O.Oe+0 
P(3.2) = O.Oe+0 
P(3.3) = t.OeO 

RETURN 

END 

* 

SUBROUTINE BNDARY(NPDE,T,U.UX,IBND, BETA, GAMMA, IRES) 

* .. Scalar Arguments .. 

DOUBLE PRECISION ALPHA 
DOUBLE PRECISION T 
INTEGER IBND, IRES, NPDE 

* .. Array Arguments .. 

DOUBLE PRECISION BETA(NPDE), GAMMA(NPDE), U(NPDE), UX(NPDE) 
DOUBLE PRECISION Cgf.cgi,Csmax,adcet,Fmass,outrad,innrad, 

+ Volflow,Eps,Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid, 

+ Deff,Vmf,Diff,bp,csv,temp,dia,Dz,yi, 

+ Gden,dens,Cpg,Cps,Ya,pi,aprime,Vel,esa, 

+ Vmw,stick,Rcons, AS,BS,CS,DS,ES,GS,FS,decet,aec,HS 

* .. Common blocks .. 

COMMON /VBLE/ALPHA 

COMMON/prop/Cgf,cgi,Csmax,Fmass,OL|trad,innrad.Volflow,Eps, 

1 Lbed,surarea,fa,Rpore,Cfs,diaf,Pvoid,aprime,Dz,yi, 

2 Deff,Vmf,Diff,bp,csv,temp,dia, VeI.esa 
COMMON/prop1/Gden,dens,Cpg,Cps 
COMMON/prop2A/mw,stick,pi,Rcons 
COMMON/prop3/decet,aec,adcet 
COMMON/Prop4/AS,BS,CS,DS,ES,GS,FS,HS 

* .. Executable Statements .. 

IF (IBND.EQ.O) THEN 
BETA(1) = 0.0 
BETA(2) = 0.0 
BETA(3) = 0.0 
GAMMA(1) = U(1)-1.0 
GAMMA(2) = U(2)-1 .0 

GAMMA(3) = U(3)-{aec*Cgf)**(1.0/fa)/((aec*CgO**{1.0/fa)+1.0) 


ELSE 



BETA(1) = 1.0e+0 
BETA(2) = 1.0e+0 
BETA(3) = 1 .Oe+0 
GAMMA(1) =0,0e+0 
GAMMA(2) =0.0e+0 
GAMMA(3) =0.0e+0 


END IF 

RETURN 

END 



